AN EXPERIMENTAL STUDY OF TAYLOR-GOERTLER 
VORTICES IN A CURVED RECTANGULAR CHANNEL 



Robert Joe McKee 



u*tr S '00! 
fo. m i 93940 



Monterey, California 





AN EXPERIMENTAL STUDY OF TAYLOR-GOERTLER 
VORTICES IN A CURVED RECTANGULAR CHANNEL 

by 

Robert Joe McKee 

Thesis Advisor: M. D. Kelleher 

June 1973 



Approved puhLLc. fizZoAML} dii&ubution wiZimi£e,d . 



An Experimental Study of Taylor-Goertler 
Vortices in a Curved Rectangular Channel 



by 

Robert Joe McKee 

Lieutenant, United States Naval Reserve 
B.S., University of California, Santa Barbara, 1968 

Submitted in partial fulfillment of the 
requirements of the degrees of 

MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

and 

MECHANICAL ENGINEER 
from the 

NAVAL POSTGRADUATE SCHOOL 
June 1973 



" 7 % S. 5 /*5 







Library 

^val Postgraduate 
Monterey, Cai-fom.. 



ABSTRACT 

A rectangular cross section channel, with both a straight 
and a curved portion, was designed and built to create the 
laminar secondary Tay lor— Goertler vortex flow. An aerosol 
was used to visualize these flow patterns, and photographs 
of the results are presented. A hot wire anemometer inves- 
tigation of one location in the curved section was conducted 
to obtain the mean velocity profiles and turbulence levels. 
The pressure drop along both the straight and the curved 
sections was measured and compared. A constant temperature 
wall heater was also designed and installed in the straight 
section of the channel. Approximate values of the heating 
losses were obtained and the overall and local heat transfer 
coefficients were calculated for the straight portion of the 
channel. It was concluded that the vortices develop with 
both velocity increases and distance along the channel, and 
that the increase in pressure drop due to the vortices 
varied from near zero to approximately 20% with increasing 
Reynolds numbers. 



\ 



2 



TABLE OF CONTENTS 

I. INTRODUCTION 11 

A. DESCRIPTION OF VORTICES 11 

B. BRIEF HISTORY 13 

II. NATURE OF THE PROBLEM 17 

A. INTENT OF THIS STUDY 17 

B. DESIGN REQUIREMENTS 18 

C. DESCRIPTION OF APPARATUS 2 4 

III. FLOW VISUALIZATION 39 

A. EXPERIMENTAL CONSIDERATIONS 39 

B. PATTERNS AND PICTURES OF THF FLOW 41 

IV. VELOCITY AND TURBULENCE MEASUREMENTS 51 

A. EXPERIMENTAL PROCEDURES 51 

B. VELOCITY PROFILES 53 

C. TURBULENCE LEVELS . 62 

V. PRESSURE DROP MEASUREMENTS 72 

A. EXPERIMENTAL PROCEDURES 72 

B. PRESSURE DROP 74 

VI. HEAT TRANSFER MFASUREMFNTS 85 

' A. EXPERIMENTAL PROCEDURES 85 

B. RESULTS 89 

VII. CONCLUSIONS 10 2 

VIII. RECOMMENDATIONS 104 

APPENDIX A: ERROR ANALYSIS 106 

APPENDIX B: COMPUTER RESULTS FOR VELOCITY 

AND TURBULENCE 10 8 



3 



APPENDIX C: CALCULATIONS AND COMPUTER RESULTS 

FOR PRESSURE DROP 119 

APPENDIX D: COMPUTER RESULTS FOR HEAT TRANSFER 129 

LIST OF REFERENCES 158 

INITIAL DISTRIBUTION LIST 161 

FORM DD 1473 162 



4 



LIST OF TABLES 



Table 



I. PRESSURE DROP PER FOOT FOR TWO FOOT 

INTERVALS CALCULATED BY LEAST SQUARES 78 

II. PRESSURE DROP PER FOOT FOR ONE FOOT 

INTERVALS CALCULATED BY LEAST SQUARES 80 

III. PRESSURE DROP PER FOOT FOR TWO FOOT 

INTERVALS MEASURED VALUES 81 



5 



LIST OF FIGURES 

Figure 

1. Sketch of Taylor-Goertler Vortices 12 

2. Reynolds Number vs. Dean Number 21 

3. Cross Section of Channel 25 

4 . Arrangement of Apparatus 27 

5. Aerosol Generating Equipment 29 

6. Lighting and Camera Arrangement 31 

7. Hot Wire Anemometer Arrangement 32 

8. Pressure Drop Measuring Equipment 34 

9. Heating Apparatus 37 

10. Photographs of Vortices Re = 609 45 

11. Photographs of Vortices Re = 875 46 

12. Photographs of Vortices Re - various 47 

13. Photographs of Vortices Re = 1033 48 

14. Photographs of Vortices Re = 1166 49 

15. Photographs of Vortices Re = 14 57 30 

16. Velocity Profiles Re = 609 55 

17. Velocity Profiles Re = 875 56 

18. Velocity Profiles Re = 1166 57 

19. Velocity Profiles Re = 1298 58 

20. Velocity Profiles Re = 1457 59 

21. Velocity Profiles Re = 1722 60 

22. Velocity Profiles Re = 1934 61 

23. Typical Position of Hot Wire * 63 



6 



24. Turbulence Intensity Re = 609 64 

25. Turbulence Intensity Re = 875 65 

26. Turbulence Intensity Re = 1166 66 

27. Turbulence Intensity Re = 1298 67 

28. Turbulence Intensity Re = 1457 69 

29. Turbulence Intensity Re = 1722 70 

30. Turbulence Intensity Re = 1934 71 

31. Pressure Below Atmospheric vs. 

Pressure Taps 75 

32. Percentage Increase in Pressure Drop 

vs. Reynolds Number 82 

33. Pressure Drop Coefficient vs. 

Pressure Taps 84 

34. Heat Transfer Coefficient Tv? = 125 95 

35. Heat Transfer Coefficient Tv; = 133 96 

36. Heat Transfer Coefficient Tw - 143 97 

37. Heat Transfer Coefficient Tw = 139 98 

38. Heat Transfer Coefficient Tw = 131 99 



7 



TABLE OF SYMBOLS 



A area 

Ac cross sectional area 

Aw area of the wall heater 

As area of a single heater strip 

a half the height of the channel 

b half the width of the channel 

Cp pressure drop coefficient 

cp specific heat 

De Dean number 

Dh hydraulic diameter 

d spacing of channel wall 

Go Goertler number 

Gr Grashoff number 

H local heat transfer coefficient 

Hm mean heat transfer coefficient 

Ho overall heat transfer coefficient 

k thermal conductivity 

L a distance along the channel in the flow direction 

Nu Nusselt number 

Pr Prandtl number 

p pressure, a relative measured value 

Qc heat energy transferred by conduction and convection 

Qr heat energy transferred by radiation 

Re Reynolds number based on hydraulic diameter 



8 



Ri 



Ro 

T 

Ta 

Tb 

T£m 

Tw 

U 

• 

V 

6 

e 

P 

a 

P 

v 



radius of inner wall or cylinder 
radius of concave wall 
temperature 
Taylor number 

bulk temperature of a fluid 
log mean temperature difference 
average temperature of the wall heater 
mean velocity of flow 

I 

volumetric flow rate 
boundary layer thickness 
emissivity 
mass density 

Stef an-Boltzmann Constant 
viscosity 

kinematic viscosity 



9 



ACKNOWLEDGEMENT 



The author wishes to express his appreciation to all 
those who made this work possible. It is impossible to give 
a complete list of people to whom the author is indebted but 
it is hoped that by naming those who deserve a special note 
of thanks, all will know that their efforts were appreciated. 

The guidance and advice of the author's thesis advisor, 
Professor M. Kelleher, and his efforts made this undertaking 
possible. The personnel of the Mechanical Engineering Shop 
and particularly Messrs. K. Mothersell, G. Bixler, and 
T. Christian deserve a special thanks for their timely and 
careful assistance in preparation of the experimental appa- 
ratus. The author is very grateful for the assistance of 
Mrs. V. Culley of the Mechanical Engineering Office in obtaining 
material and completing administrative requirements. Final- 
ly, to my wife Catherine, for her patience through repeated 
manuscript typings and for her support throughout the study, a 
heartfelt thank you. 



10 



INTRODUCTION 



I . 

A. DESCRIPTION OF VORTICES 

The phenomenon best known as Taylor-Goertler vortices is 
caused by a body force instability in fluid flows. Speci- 
fically, Taylor-Goertler vortices are a secondary flow 
induced by centrifugal forces. In a channel, such as used 
in this investigation, the fluid near the center of the 
channel is subjected to larger centrifugal forces than the 
slower moving fluid near the concave wall. Thus, the ten- 
dency is for the fluid in the center of the channel to move 
outward toward the wall. The fluid near the wall is unable 
to resist this action and must move in the spanwise direction 
and then radially inward replacing the central fluid. Once 
in the center of the channel, this fluid obtains the higher 
flow velocity, and the tendency to move toward the concave 
wall. This cyclic motion forms the rotating Taylor-Goertler 
vortices . 

Such vortices, caused by centrifugal forces, occur in 
many cases of curved flow. Taylor-Goertler vortices are 
laminar vortices with their axes in the direction of the 
main flow, and with secondary velocities in both perpendic- 
ular directions. Such spiral vortices occur in counter 
rotating pairs. There is frequently a regularly spaced 
cellular structure associated with these vortices. A sketch 
of this type of flow pattern is shown in Figure 1. 
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Figure 1. 



Sketch of Taylor-Goertler Vortices. 
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B. BRIEF HISTORY 



The first person to consider the instability of curved 
flows was Lord Rayleigh who found a stability criterion for 
inviscid fluid rotating symmetrically about an axis [Ref. 1] . 
His criterion for a stable flow required that the product of 
the local circumferential velocity and the local radius of 
curvature either increase or at least remain constant as the 
radius increases. An example of an unstable flow described 
by Lord Rayleigh is the flow between a fixed outer cylinder 
and a concentrically rotating inner cylinder. In 1923, 

G. I. Taylor published an extensive analytic and experimental 
study of such unstable Couette flows of viscous fluids 
between rotating cylinders [Ref. 2]. The Taylor number, 
defined as 



was shown to be a characteristic parameter of such flows with 
a critical value for the onset of the instability. The 
secondary flow observed between the cylinders when the outer 
was fixed and the inner rotated has been shown to be vortices 
whose axes were in the circumferential direction. These 
vortices were named as Taylor vortices. In recent years, 
there have been a large number of investigations’ of cylindri- 
cal Couette flow. The torque required to rotate a cylinder 
in the presence of Taylor vortices was measured by P. Castle, 
et al. [Ref. 3]. The same type of instability in the boundary 
layer along concave walls was first investigated by 
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H. Goertler in 1940 [Ref. 4]. The Goertler number, which 
is similar to the Taylor number, is a parameter characteris- 
tic of curved boundary layer flows and is defined as 




where 6 is the boundary layer thickness . When the Goertler 
number exceeds a critical value, the presence of Goertler 
vortices is observed. H. W. Liepmann [Ref. 5] was among the 
investigators to continue and verify Goertler' s experimental 
work. The approximate analytic results obtained by Goertler 
have been verified with an exact solution by the German 
investigator G. Hammerlin as reported by H. Schlichting 
[Ref. 6]. These results were also verified with an extensive 
numerical solution completed by A. M. 0. Smith [Ref. 7]. 

The problem of instability of viscous fluids flowing in 
a curved channel was first considered by W. R. Dean in 1928 
[Ref. 8]. Dean found a parameter, similar to the Taylor and 
Goertler numbers , which governs the onset of Tay lor-Goertler 
vortices in curved channels. The analytic work done by Dean 
has been verified by W. H. Reid [Ref. 9] among others. The 
experimental work on secondary flows in channels seems very 
limited. Until just recently, the possible effects of 
Tay lor-Goertler vortices on heat transfer had not been 
studied. In 1965, L. Persen [Ref. 10] considered the effect 
of Goertler type vortices on the heat transfer from a wall 
for the special cases of both a very high and a very low 
Prandtl number. The first experimental work with the effect 
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of Tay lor-Goertler vortices on heat transfer in a boundary 
layer was reported by P. McCormack et al. [Ref. 11] in their 
1970 study. With some limiting assumptions, a numerical 
solution for the heat transfer in curved rectangular chan- 
nels was reported by K. Cheng and M. Akiyama in 1969 
[Ref. 12], Analytic and experimental results for a fully 
developed, constant wall heat flux, square cross section 
curved channel were obtained in 1970 by Y. Mori, et al . 

[Ref. 13] . 

Although Tay lor-Goertler vortices are a laminar flow, 
they have secondary velocity components and are more complex 
than one would expect. Furthermore, as indicated by H. W. 
Liepmenn [Ref. 5], these vortex flows are potentially a key 
to understanding transition to turbulent flow. Tay lor- 
Goertler vortices have similarities to other vortex flows 
such as wing tip vortices and the vortex rolls in forced 
convection heating of fluid layers described in the paper by 
M. Akiyama et al. [Ref. 14]. The striations seen at stag- 
nation points on bluff bodies and the cross hatching on 
reentry bodies have also been attributed, at least in part, 
to the presence of Tay lor-Goertler vortices, as reported by 
L. Persen [Ref. 15] . 

There are many possible applications that could result 
from investigation and understanding of these vortex phenom- 
enon. Potential application include improved cooling of 
turbine blades and other external surfaces, reduction of 
pressure losses in bends, and improved mixing in laminar 
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flows. Heat exchangers which could utilize the low pressure 
drop in laminar flow as well as the enhanced heat transfer 
in the presence of vortices would be a significant advance- 
ment in heat transfer. 
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II. NATURE OF THE PROBLEM 



A. INTENT OF THIS STUDY 

Because future application of the Tay lor-Goertler vortex 
phenomenon may take place in closed channels, it was intended 
to study the effect of these secondary flows in a curved 
channel of rectangular cross sectional area. To investigate 
Tay lor-Goertler vortices in a curved rectangular channel, the 
presence of the vortices had to be verified. It was intended 
to visualize the flow in order to confirm the presence of 
the vortices. How rapidly the vortices would form and to 
what extent they would fill the channel was not known. It 
was intended to use the flow visualization technique to 
study, as far as possible, the size, the form, and the 
development of the vortices. 

Other planned measurements included a hot wire anemometer 
check of the turbulence level to insure that the flow was 
laminar. The velocity profiles in the channel were also to 
be obtained to determine what effect the secondary flow would 
have on the mean velocity. The additional velocity com- 
ponents in Tay lor-Goertler vortices produce Reynolds stress 
not present in other laminar flows. It was planned to 
measure the pressure drop along both the straight and 
curved portions of the channel used in this investigation 
to determine the effect, if any, of these Reynolds stresses 
on the pressure drop. An increase in the pressure drop 
would be anticipated because it is known that additional 
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energy is required to maintain the secondary flow. This 
energy is supplied to the Taylor-Goertler vortices by an 
interaction of the mean and secondary velocity components. 
The works of G. I. Taylor [Ref. 2] and P. Castle et al. 

[Ref. 3] have shown the increase in torque required to 
rotate a cylinder in the presence of Taylor vortices. 

The Taylor-Goertler vortices would also be expected to 
effect the heat transfer from the concave wall. This effect 
would result from the secondary velocity components trans- 
porting heated fluid from the wall towards the center of the 
channel, and carrying cooler fluid from the center to the 
wall. It was intended to study the requirements for measur- 
ing the effects of this action on the heat transfer from 
the concave wall. 

B. DESIGN REQUIREMENTS 

In order to study Taylor-Goertler vortices, an apparatus 
had to be built which would give the values of the non- 
dimensional parameters required to produce the secondary 
flow. The limit of stability of the basic laminar flow in 
a curved channel against the formation of Taylor-Goertler 
vortices is well known, and is presented as a critical Dean 
number by W. Reid [Ref. 9]. If the Dean number, defined as 



exceeds a value of 36, it is presumed that Taylor-Goertler 
vortices will be present. Furthermore, it is known that the 
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rate of development of the vortices will be somewhat depen- 
dent on how much greater than 36 the Dean number becomes. 

It is known that the size of the vortices will depend in 
some manner on the dimensions of the channel. Investigators 
of Tay lor-Goertler vortices on concave walls have reported 
vortices from the size of the boundary layer to ten times 
larger, dependent on geometric factors [Ref. 7]. The Dean 
number had to be greater than the critical value of 36 and 
be adjusted within a workable range. One hundred fifty was 
taken as the upper limit of this range, as it was assumed 
the vortices would be fully developed by at least this value. 

The minimum critical Reynolds number, based on hydraulic 
diameter, at which turbulent flow is observed in a parallel- 
plate channel is given by G. Beavers, et al. [Ref. 16] as 
Re = 2200. Several investigators including H. Liepmann 
[Ref. 5] and R. Nunge [Ref. 17] have reported that the 
decrease in the critical Reynolds number due to curvature is 
small, especially for a large radius of curvature. There- 
fore, the maximum Reynolds number, based on hydraulic dia- 
meter, selected for this investigation was Re = 2000. The 
Reynolds number had to be high enough for the flow velocities 
to be measured accurately. A minimum of Re = 100 was used 
for the design of the experimental apparatus. 

The design procedure involved checking the effects of 
several factors on the range of the nondimensional variables. 
The working fluid, the channel height, and the flow rate 
selected, each have an effect on both the Reynolds number 
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and the Dean number. The radius of curvature also effects 



the Dean number. Trial and error was used to select the 
height and width of the channel's cross section, and the 
following inequality in the Reynolds number was checked to 
give the range of velocities 

100 < Re < 2000 

The inequality in the Dean number given below was then 
checked to give the range of the radius of curvature. 

36 < De < 150 

When the velocity range was practical to obtain and easy to 
measure, and when the radius of curvature was reasonable for 
construction, the design variables were considered feasible. 
The resulting relationship between the Reynolds number, 
based on hydraulic diameter, and the Dean number with the 
ratio of radius of curvature to plate spacing as a parameter 
is shown in Figure 2. The selected design is represented by 
the line R/d = 48. 

The design of the apparatus for this investigation was 
also constrained by time and cost limitations. The material 
and equipment to be used had to be workable, available with- 
out long delay, and inexpensive. Water and air were the 
only working fluids considered as practical. The flow 
velocities at which water satisfied the Reynolds number and 
the Dean number requirements were too low to measure by 
ordinary means. Furthermore, water requires more elaborate 
construction to prevent leakage from the channel. For these 
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Dean number, De 



reasons, air was selected as the working fluid. For design 
purposes, the value of the fluid properties at estimated 
operating temperatures were taken from reference tables. 

The length of the entrance region in the design channel was 
checked to see that it would not be excessive. To insure 
smooth entering flow, the entrance contraction section was 
designed in accordance with considerations given by M. Cohen 
and N. Ritchie [Ref. 18]. The intent to visualize the flow 
required the use of a workable, transparent, inexpensive 
material for the walls of the channel. Plexiglas was 
selected for this purpose. 

It was hoped that the results obtained for both pressure 
drop and heat transfer would be comparable to infinite 
parallel plate solutions, so a large aspect ratio was chosen. 
The pressure drop for flow between parallel plates is known 
from an exact analytic solution and is not of a large magni- 
tude. This magnitude was compared with the sensitivity of 
the laboratory micromanometer available. It was decided 
that the pressure drop could be measured, but only if a 
sufficiently long test section was built. No exact solution 
exists for the pressure drop in a curved channel in the 
presence of Tay lor-Goertler vortices. However, the studies 
by R. Nunge [Ref. 17] and by K. Cheng and M. Akiyama [Ref. 12] 
indicates a rather small increase in the friction factor in 
curved channels of large aspect ratio. It was assumed that 
the curved test section would be sufficiently long to allow 
measurement of the pressure drop. 
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To measure the effects of Tay lor-Goertler vortices on 
the heat transfer from a constant temperature concave wall 
of a curved channel requires comparative data for one wall 
of a flat plate channel. Heat transfer data for a flat 
plate channel with one constant temperature wall is very 
limited. Limiting values of the Nusselt number and other 
comparative values can be obtained from R. Shah and A. Lon- 
don's report [Ref. 19]. The design of a heater plate to 
measure this type of heat transfer presented special require- 
ments. As with any flat plate heater, more power would be 
convected away from the upstream edge than from points 
farther downstream. Therefore, to obtain an uniform wall 
temperature, by electrical heating, separate resistance 
elements with individual voltage controls were required. A 
multiple channel regulated power controller was not immedi- 
ately available, but it was determined that such an instru- 
ment could be constructed. Grade A nichrome was selected 
for the heating elements because of this metal's large 
specific resistance, small temperature coefficient of resis- 
tance, and relative low cost. The many power leads and 
temperature recording leads would have to reach the heating 
element through the channel wall. For this reason a castable 
plastic resin was selected for constructing the heater sec- 
tion. Although only a flat heater plate was to be constructed 
for the first tests, design considerations were given to 
adapting this construction method to a heater plate for the 
curved section. Other aspects of the heated section con- 
struction will be discussed later. To insure that the heat 
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transfer data v/ould be for forced convection, and would not 
contain the effects of free convection, the value of the 
governing dimensionless group was checked against the 
requirement 

Gr * Pr * — < 1000 
L - 

given by J. Holman [Ref. 20]. It was found that this 
requirement imposed a significant temperature difference 
limitation for water as the working fluid, but allowed a 
large temperature range for the selected working fluid, air. 

The effects of the many design requirements discussed 
above was to make the design problem somewhat difficult. 
However, the v/orking fluid, the geometry, and the flow 
velocities to be studied are of values frequently encountered 
in practice. 

C. DESCRIPTION OF APPARATUS 

The cross section of the channel was 0.250 inches high 
and 10.0 inches wide. It had an aspect ratio of 40 and a 
cross sectional area of 0.017361 square feet. The hydraulic 
diameter of this channel v/as 0.040650 feet. The channel was 
made of two one-quarter inch thick sheets of plexiglas 
separated by rolled metal spacers which also formed the 
sides of the channel as shown in Figure 3. Steel was used 
for the spacers in the straight section of the channel and 
for added flexibility, aluminum was used in the curved por- 
tion. A straight portion of this channel, consisting of a 
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two foot long entrance section and a two foot long test 
section, preceded the curved section. The radius of curva- 
ture of the interior concave wall of the curved section was 
one foot. The curved section was formed by heating the 
plexiglas and bending it over a carefully constructed wood 
frame. This wood frame was also used to support the entire 
channel. The curved portion of the channel completed a 
180 degree turn and was followed by a short straight section. 

The flow of air entered the channel through a contrac- 
tion section attached to the beginning of the entrance region 
and covered by a cheesecloth screen for most of the experi- 
mental work. The flow left the channel through an exhaust 
nozzel which was connected to flexible tubing of one inch 
inside diameter . The tubing led the flow through a Fischer 
and Porter Company variable area flow meter, model number 
10 A3565. This rotometer had a 100% full scale flow rate 
of 11.1 standard cubic feet of air per minute and an accuracy 
of ±0.5% of full scale. The flow was drawn through the flow 
meter, the flexible tubing, and the channel by an electri- 
cally driven Cadillac, model G12 , centrifugal blower. The 
blower speed was controlled by varying motor voltage with a 
variac. The voltage was supplied by a Sorenson, R 1050, 

A. C . voltage regulator. The general arrangement of this 
apparatus is shown in Figure 4. 

The idea of using an aerosol for visualizing the flow 
was taken from a paper by O. Griffin and C. Votaw [Ref. 21]. 
The advantages and ease of using an aerosol are derived from 
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Figure 4. Arrangement of Apparatus. 



the non-toxic and non-corrosive properties of the DOP , 
di (2-ethylhexy 1) -phthalate , chemical used. The equipment 
used to make the aerosol was built in accordance with a 
United States Naval Research Laboratory report by W. Echols 
and J. Young [Ref. 22] and personal communications with 
0. M. Griffin. Air from a regulated and filtered supply was 
piped to an atomizer nozzle submerged in a jar of DOP. In 
addition to the pressure regulator, there was a throttling 
valve on the air line and the depth of submergence of the 
nozzle could be adjusted. The aerosol from the top of this 
jar was collected and led to a jet impactor in a second jar. 
The effect of the jet impactor was to remove excessively 
large particles and to make a smaller, more uniform aerosol 
particle size. The aerosol was then led to a settling 
chamber where the pressure was reduced to nearly atmospheric 
level by bleeding off the excess aerosol. Flexible injec- 
tion tubing of 0.250 inch inside diameter connected the 
settling chamber to the contraction section where the 
aerosol entered the channel. The aerosol generating equip- 
ment is pictured in Figure 5. The character of the aerosol 
used was such that visualization of the flow was facilitated 
by illuminating slits along the channel and viewing them at 
approximately 90° to the light. A single 650 watt Colortran 
light, model 100-071, was placed on the side of the channel 
opposite the viewer. The light was collimated with a paper 
cone lined with aluminum foil. Flat black paper was taped 
to the outside of the channel walls to make four narrow 
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Figure 5. Aerosol Generating Equipment. 
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slits. The light beam was directed through one of these 
slits so that a small area of the flov; could be visualized. 

A Nikomat 35 mm camera and a Nikon 55 mm f3.5 Macro-Nikkor 
lens was used to obtain photographs of the flow patterns. 

A sketch of the lighting and camera arrangement is shown in 
Figure 6. The visualization techniques and procedures will 
be discussed more fully when the photographs are presented. 

In addition to visualizing the flow it was planned to 
measure the velocity profile. Due to the size of the chan- 
nel, a sub miniature hot wire anemometer probe, Thermo- 
Systems Inc. model 1279, which had a 90° bend in the supports 
was used to measure the turbulence levels . This probe was 
located three inches from the center line of the channel and 
20 inches of arc length downstream from the start of the 
curved portion of the channel. The probe was placed in a 
plug and inserted through the convex wall of the channel. 

The inside surface of this plug fit flush with .the wall of 
the channel, and had the same radius of curvature. The 
0.035 inch diameter body of the hot wire anemometer protruded 
from the plug and was fastened to a micrometer barrel which 
was mounted on a machined block attached to the outside of 
the channel. In this way, the movement of the wire from 
very close to the concave wall of the channel to near the 
convex wall could be controlled and measured. This arrange- 
ment is shown in Figure 7. The leads from this probe were 
connected to a Thermo-System Inc. hot wire anemometer 
system, model 1050, including a Hewlett-Packard RMS Meter. 
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igure 6. Lighting and Camera Arrangement. 




Figure 7. Hot Wire Anemometer Arrangement. 



32 



A Thermo-Systems Inc. calibrator, model 1125, was used to 
calibrate the hot wire anemometer. 

In order to measure the pressure drop along the length 
of the channel, pressure taps and a micromanometer were 
required. Pressure taps were located every three inches 
along the entire length of the channel, starting two inches 
from the beginning of the entrance section and alternately 
set 0.250 inches to the right or left of the center line. 
Along the curved section, the three inch spacing was 
measured on an arc of one foot radius, resulting in a total 
of 28 pressure taps with number 17 being two inches into 
the curved section. These pressure taps consisted of a 
0.040 inch diameter hole drilled into the channel with a 
short plexiglas tube glued into a 0.250 inch hole counter- 
bored over the smaller holes. Flexible tubing with pressure 
fittings was connected to these short tubes to facilitate 
connecting and disconnecting the micromanometer.. The instru- 
ment used in this investigation was an E. Vernon Hill and 
Company Type C micromanometer with a range up to two inches 
of water and a sensitivity of 0.001 inches of water. This 
micromanometer is of the re-zeroing type where the reservoir 
is raised to return the meniscus to the same level it had 
before the pressure was applied, and the amount the reser- 
voir is raised indicates the pressure. A photograph of the 
channel and the pressure measuring apparatus is shown in 
Figure 8. 
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The bulk temperature of the air was measured six inches 
after the entrance of the channel, and three inches after 
the heating plate, by four thermocouples spaced across the 
channel at each location. Copper-constantan thremocouples 
were used. They were located at different heights in the 
channel and connected in parallel to give the average bulk 
temperature. An ice bath reference junction was provided, 
and these thermocouples and the others used in this experi- 
ment were connected through a switch to a Dymec digital 
voltmeter, model 2401B. These thermocouples and a sample 
of the others used in this study were calibrated in a 
Rosemont Laboratory Calibration System. In this constant 
temperature bath, the temperature was determined by a 
platinum resistance thermometer. 

A wall heater was constructed for the straight test 
section of the channel. This heater was formed of forty 
0.004 inch thick, 0.125 inch wide strips of nichrome run 
across its 11.5 inch width. These metal ribbons were 
placed in a mold by hand and spaced approximately 0.010 
inches apart. The streamwise length of the plate was 
5.50 inches. A copper-constantan thermocouple was welded 
to the center of each of the nichrome strips. Power leads 
were connected to one side of each strip four inches from 
the center of the plate. On the other side of the center- 
line, the strips were connected together in pairs to make 
20 electric circuits. The 20 heaters were individually 
supplied by a specially constructed, regulated 
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high current D.C. power supply and a 20 channel power con- 
troller. The D.C. power supply had a 200 amp capacity and 
a maximum RMS ripple of one millivolt. The power controller 
was a 20 channel potentiometer biased series regulated 
voltage source with a temperature compensated reference 
voltage. The components of this power controller were 
mounted on a common heat sink in order to improve the 
stability of operation. The power supply and controller 
are shown in operation in Figure 9. When all the power 
leads and thermocouple leads were in place, APCO 210 epoxy 
resin with APCO 180 hardener was poured into the mold. This 
thermal set plastic was cured at room temperature for 
approximately 24 hours. When the heater plate was removed 
from the mold, it was found to be warped. A metal frame was 
constructed, securely attached to the outer edges of the 
0.250 inch thick casting and adjusted to produce a completely 
flat plate. A further problem experienced with the casting 
involved some of the nichrome strips near the trailing edge 
of the plate that had pulled away from or sunk into the 
plastic during drying. This resulted in rough spots on the 
surface of the plate. The loose strips were cut and removed 
from the heater, and the low spots, from the 33rd strip to 
the end of the plate, were filled with epoxy. The plate was 
then lightly sanded, resulting in a smooth surface. The 
first 16 heating circuits remained useable, which resulted 
in a heated area 8.125 inches wide by 4.384 inches in the 
flow direction. 
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After the other experimental work was completed, a 
5.50 inch piece of the upper wall of the channel, between 
the two foot ten inch station and the three foot 3.50 inch 
station, was removed. The heater was placed into this cut, 
becoming part of the straight wall of the channel with the 
nichrome ribbons exposed to the flow of the air. The top 
of the heater and the surrounding channel walls were then 
covered with a 0.125 inch thick layer of Johns-Manville 
Company Min-K insulation. The Min-K was covered with a one 
inch layer of glass wool. Each piece of the apparatus wa s 
cleaned and carefully prepared before the experimental work 
began. 
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III. 



FLOW VISUALIZATION 



A. EXPERIMENTAL CONSIDERATIONS 

The aerosol generating equipment was found to be very 
versatile as to the amount and quality of the aerosol pro- 
duced. After much trial and error, good quality aerosol was 
found to be that of relatively small particle size and low 
exhaust pressure. This aerosol appeared much like cigarette 
smoke when discharged into the room. A supply air pressure 
of 30 pounds per square inch, or slightly less, an atomizer 
nozzle submergence of approximately 0.50 inches, and a low 
settling chamber pressure were the conditions which produced 
this aerosol. Producing this quality aerosol and getting 
the right amount of it into the channel proved to be impor- 
tant in visualizing the flow. Viewing the aerosol at the 
proper angle relative to the incident light was also impor- 
tant. For lighting, the best arrangement was to place a 
single, strong light near the center of curvature of the 
channel and to shine it radially outward. The collimated 
beam of light was directed through narrow slits formed from 
black paper taped to the channel. If this small area of 
the flow was viewed from an angle tangent to the channel, 
the vortex patterns could be readily seen. This viewing 
location was the same as the camera location shown in 
Figure 6. Four such narrow slits were made along the curved 
section of the channel. The first slit was located 7.50 
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inches from the beginning of the curved portion of the 
channel. The second was 12.00 inches downstream from the 
first, and the third was 11.50 inches farther downstream. 

The last slit was located at the end of the curved section, 
8.0 inches from the third slit. At the location of this 
fourth slit, the channel enters the final straight section. 

To obtain photographs of the flow patterns in each of 
these slits, the camera with the 55mm macro lens was mounted 
on a tripod and moved as close to the channel as possible. 
The camera was aligned with the axes of the vortices by 
placing the camera on a tangent to the channel. In the case 
of the first and second slits, the flow was coming toward 
the camera, and in the case of the third and fourth slits, 
the flow was away from the camera. Due to the many reflec- 
tions observed from the surfaces of the plexiglas, the room 
was darkened while the photographs were taken. A Gossen 
"Luna-Pro" exposure meter with a variable angle spot meter 
attachment was used to assist in obtaining the proper 
exposure settings. Under these difficult conditions, several 
attempts were required before satisfactory photographs were 
obtained. 

Other arrangements for viewing the flow were also tried. 
Looking radially inward through the channel toward a less 
intense light located near the center of curvature showed a 
different type of pattern. The patterns seen were streaks 
or lines in the direction of the flow. These patterns were 
similar to those seen by other investigators including 
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L. Persen [Ref. 15] . Several methods of injecting the aero- 
sol into the f lev were also tried. The methods involving 
small tubes or fittings did not work because sufficient 
aerosol would not pass through these passages. The results 
obtained from the other methods were for the most part the 
same as for the method chosen. In a few cases, the early 
development patterns appeared differently, depending on how 
the aerosol entered the channel. The fully developed vortex 
pattern, however, always appeared the same. 

B. PATTERNS AND PICTURES OF THE FLOW 

The pictures of the vortex patterns are presented in 
Figures ten through 15. Belov; each picture is a slit number, 
a Reynolds number and a Dean number indicating the location 
of the slits, the flow rate and the value of the character- 
istic parameter for secondary flow respectively. For the 
channel used in this investigation, the Dean number is 
related to the Reynolds number by 

De = 0.073973 * Re. 

The line below the vortex patterns in these enlarged pictures 
is the concave wall of the channel. The line of reflected 
light above the vortices is the inner wall of the channel 
which is 0.25 inches from the concave wall. The photographs 
taken at slit four appear slightly different than the others 
because the lighting angle and therefore the reflections 
from the plexiglas were slightly different. The pictures in 
each figure are arranged with slit one at the top of the page 
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and the downstream slits in order below the first. The 
Reynolds numbers of the pictures in each figure are intended 
to be of the same magnitude. There is not complete repre- 
sentation of the same Reynolds numbers or an even sampling 
of the slit location because of the moderately high failure 
rate of obtaining pictures. In spite of this difficulty, 
the photographs presented are indicative of the Taylor- 
Goertler vortices as observed. The figures are in order of 
increasing average Reynolds number. 

A number of distinct patterns appear repeatedly in these 
pictures. At the lowest velocities, the aerosol remains 
near the concave wall in a pattern that is not clearly a 
fully developed vortex flow. Examination of these patterns 
in Figures 10a, 10b, 11a, and 12a show that the ends are 
curled up. This same pattern has been observed farther 
downstream in slit three at flow rates lower than pictured 
here. This pattern occurs predominately near the beginning 
of the curved section. From extended observation of both 
this pattern and the motion at the curled end, it is felt 
that they represent the beginning of a developing weak vortex. 

The next form considered looks much like a "mushroom." 

This pattern is exemplified by Figures 12b, lib and the right 
side of 10c. The appearance of two different patterns under 
the same conditions, as in Figure 10c, was observed occasion- 
ally. Such occurrances are possibly explained in terms of 
the actual finite width of the channel which produces a 
stronger velocity gradient near the side walls than in the 
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middle. A difference in the placement of the aerosol injec- 
tion tubes could also contribute to the appearance of dif- 
ferent patterns at the same flow rate. Furthermore, as the 
vortices develop and occupy more space in the channel, some 
will be distorted to make room for the others. At higher 
Reynolds numbers, the top of the curl pattern, as in 
Figures 13a, 14a, and 15a, could be considered a "mushroom" 
on a small scale. The "mushroom" shape occurs farther down- 
stream or at higher velocities than the first pattern 
discussed. This shape always proceeds in distance or flow 
rate, the first clear vortex patterns. The "mushroom" has 
been observed to lead directly into distinct vortices and 
therefore should be considered as a developing vortex. 

The next pattern is the round form shown in Figures 11c, 
13b, 14b, and 15b. These forms represent a counter rotating 
pair of vortices with flow leaving the concave wall in the 
center and curling back toward the wall on either side. 

These round vortex patterns do not appear to occupy all the 
space in the channel. In Figure 13b there appears to be 
space above the vortices, and in Figure 15b the space appears 
below the vortices. However, the aerosol which is what was 
photographed does not necessarily show the entire vortex 
pattern. The aerosol shows a streakline line pattern which 
may be somewhat dependent on how the aerosol enters the 
developing vortices. As the flow rate and the downstream 
distance increases, the fully developed vortex patterns 
appear. These vortices, as shown in Figures lOd, lid, and 
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12c are the full height of the channel and occupy a nearly 
cellular area. Such patterns could be considered as fully 
developed Taylor-Goertler vortices. 

As the velocities continued to increase, the vortices 
grew in width until they began to crowd one another, as seen 
in Figures 13c, 14c, 14d and 15d. If the distance or the 
Reynolds number were increased farther, an unsteady oscil- 
latory motion developed in the spanwise direction. This 
action which could be observed by eye, was difficult to show 
in a still photograph. Figure 13d is the result of such an 
attempt. The observations indicated that the sideways oscil- 
lations of the vortices may involve the absorption of one 
vortex by an adjacent growing one. It was apparent that an 
increase of oscillatory motion would have led to turbulence. 

It was evident from the photographic study that the 
vortices did require some distance to develop. Several 
different flow patterns were observed during the development. 
At high-flow rates, less distance was required for the 
vortices to become full size. There were several complex- 
ities observed in the flow patterns which could lead to 
further studies. 
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(b) . Slit 3 Re = 609 De = 45.0 



(a). Slit 1 Re = 609 De = 45.0 



(d) . Slit 4 Re = 742 De = 
Figure 10 
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(d) . Slit 4 Re = 874 De = 64.7 
Figure 11 
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(c) . Slit 3 Re = 927 De = 68.6 




(d) . Slit 4 Re 



874 De = 64.7 



Figure 12 



47 







(b) Slit 2 Re = 1033 De = 76.4 



(c) Slit 3 Re = 1033 De = 76.4 



(d) Slit 4 Re = 1033 De = 76.4 
Figure 13 



(a). Slit 1 Re = 1033 De = 76.4 
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Figure 14 
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(a) Slit 1 Re = 1457 De = 107.8 





(b) Slit 2 Re = 1245 De - 92.1 



(c) Slit 2 Re = 1457 De = 107.8 




(d) Slit 3 Re = 1298 De = 96.0 



Figure 15 
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IV. VELOCITY AND TURBULENCE MEASUREMENTS 



A. EXPERIMENTAL PROCEDURES 

The sub miniature hot wire anemometer probe to be used 
to obtain the mean velocity profiles was connected to the 
anemometer. The cold resistance of the wire was measured 
as 7.90 ohms. Based on an overheat ratio of 1.5, the 
anemometer's operating resistance was set at a value of 
11.85 ohms. The anemometer's signal was fed into an oscil- 
liscope so that the patterns could be observed. Stability 
and trim controls were adjusted for maximum frequency 
response. The probe was centered in the chamber of the 
calibrator. A filtered air supply was connected to the 
calibrator, and the pressure drop across the calibrator's 
flow nozzle was measured by an incline manometer. With the 
anemometer switched to run, the flow rate through the 
calibrator was adjusted to different levels. The pressure 
drop in inches of water, the bridge output voltage, and the 
RMS voltage were recorded. This data was reduced by the 
calibration program shown in Appendix B and the square of 
the output voltage was plotted as a straight line against 
the square root of the flow velocity in feet per second. 

The slope and intercept of that straight line were the 
values required to reduce further data, and were found to 
be 0.710 and 2.020 respectively. 
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The hot wire anemometer was removed from the calibrator 
and inserted as far into the plexiglas mounting plug as 
possible. The mounting block for the micrometer barrel was 
attached to the stem of the probe, and the distance from 
the base of the block to both the top of the plug and to 
the hot wire were measured. The difference in these machine 
vernier readings was 0.050 inches. This was as close as 
the hot wire anemometer could come to the convex wall of the 
channel. The micrometer barrel and a straight piece of 
steel wire had been used to measure the distance across the 
channel at the location of the plug. This distance was 
found to be 0.240 inches, indicating that the probe would 
have a travel of 0.190 inches across the channel. When the 
plug was inserted into the channel and the micrometer was 
attached to the mounting block, the micrometer reading was 
0.460 inches. The hot wire anemometer was carefully moved 
until the supports for the wire contacted the concave wall 
of the channel. The micrometer reading at this point was 
0.270 inches, verifying the 0.190 inches of travel. The 
wire itself was very close to, but not touching, the wall. 
The data was taken between these micrometer settings of 
0.270 inches and 0.460 inches. 

Before the data was taken, the channel, the blower, and 
the electronic instruments were allowed a warm up period of 
about one hour. After the warm up, the flow rate was 
adjusted and a minimum of 30 minutes was allowed for steady 
state to be achieved. At each flow rate, the temperature of 
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the incoming fluid was measured and the micrometer was moved 
0.010 inch between data positions. At each position, the 
flow rate was checked to see that it was stable, and the 
micrometer setting, the output voltage, and the RMS voltage 
were recorded. This data was reduced by a computer program 
which plotted the mean velocity and the turbulence intensity 
against the distance from the concave wall. This program 
and the printed results are contained in Appendix B. 

B. VELOCITY PROFILES 

The velocity recorded and shown in the figures which 
follow was not just the mean velocity component down the 
channel, but the vector sum of the streamwise and the radial 
velocity components. The portion of the response that was 
due to the radially inward or outward flow could not be 
determined by this one set of measurements. However, at the 
flow rates of this experiment, a nearly parabolic velocity 
profile would be expected if the secondary flow was not 
present [Ref. 17] . This has also been shown for curved 
channels of square cross section [Ref. 13] . The position 
of the hot wire relative to the vortices was not determined 
during the tests. Considering the nature of vortex motion, 
a shift in the position of the vortex relative to the probe 
could have made a difference in the velocity indicated. The 
location of the hot wire anemometer was very close to the 
second slit used in flow visualization. Thus, for each flow 
rate, the vortices may be considered to have developed to a 
corresponding level. 
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Figures 16 through 22 show the mean velocity in feet per 
second versus the distance from the concave wall in hun- 
dredths of an inch. There was a maximum uncertainty of 
approximately 0.010 inches on the location of the wall. The 
uncertainty in the velocity measurement using the T.S.I. 
system was at most 0.05 feet per second. Figure 16 shows 
a parabolic shape with a small increase in velocity indicated 
near 0.060 inches. This point may be associated with an 
active point of a developing vortex. The step in the velo- 
city shown near the wall could be associated with the probe 
actually contacting the wall, or with the error in position- 
ing. In this case, and in all other graphs, the curve should 
be extended to zero at the wall. Although a vortex was most 
probably present, Figure 17 shows no irregularities in the 
velocity profile. Perhaps the probe was located near the 
center of a vortex at this time. The same comment could 
be true of Figure 18. The magnitude of the velocity is 
increasing as would be expected. 

Figure 19 shows a variation from the parabolic profile. 

As indicated by Figure 15b, the vortices were growing at 
this Reynolds number and a different portion of the vortex 
may have been acting on the wire. In Figure 20 the velocity 
profile appears much different. In this case, the radial 
velocity component was larger near the wall than in the 
center of the channel. Figure 21 is much the same as 
Figure 20. A different profile appears in Figure 22. This 
last profile could be explained in terms of the hot wire 



54 



o 



o 



O 

LD 




O 

CD 

If) 

u 

CD 

OU 

-P 

Q) 

CD 

4 -< 



O >i 

• V 
m -H 
O 
O 

r— I 

CD 

> 



CN 



I L_ 

O 

m cn 

CN 

• O 

o 



I I 

• • 

o o 

(seqouT) eouBqsxp 



j 

LD 

O 

O 



o 

o 



o 

o 



55 



Figure 16. Velocity Profile Re = 609 
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Figure 18. Velocity Profile Re = 1166. 
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Figure 20. Velocity Profile Re = 1457 
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velocity (feet per sec) 



moving through the edge of a vortex. Near the wall the 
outward radial velocity would increase the mean velocity 
measured. In the center of the vortex there may be little 
effect on the hot wire anemometer. Near the top of the 
vortex the radial component of the velocity would again 
affect the probe. At the other edge of this vortex, the 
inward movement of the fluid v/ould produce the same velocity 
profile. The length of the hot wire, 0.063 inches, would be 
sufficiently small to traverse the vortices in this manner. 

A schematic drawing of the probe in a typical position 
relative to the vortices is shown in Figure 23. 

C . TURBULENCE LEVELS 

The turbulence intensity was defined as the root-mean- 
square of the velocity fluctuation divided by the mean 
velocity. These values were calculated by the data reduc- 
tion program from the RMS meter readings and the bridge 
output voltages. The results were plotted against the 
distance from the • concave wall in Figures 24 through 30. 

The turbulence level was quite low in the center of the 
channel and higher near the walls as a result of the low 
mean velocity near the walls. 

The change in the turbulence level near the wall in 
Figure 24 was due to the step change in the velocity profile 
for that location and Reynolds number. Otherwise, Figures 
24, 25, and 26 are self explanatory. In Figure 27, the 
curve representing the turbulence intensity is not as smooth 
as the other curves. Large increases in the fluctuations of 
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Figure 23 . Typical Position of Hot Wire 
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Figure 24. Turbulence Intensity Re = 609. 
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Figure 25. Turbulence Intensity Re = 875. 
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Figure 26. Turbulence Intensity Re = 1166. 
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Figure 27. Turbulence Intensity Re = 1298 



the RMS meter started at this flow rate, resulting in an 
increased uncertainty in the level of turbulence. With the 
RMS meter reading at the 0.00350 volt level, the value would 
suddenly increase as much as 0.010 and drop back again. 

The frequency of this random increase was of the order of 
30 cycles per minute. One possible cause of these large 
jumps in the local velocity would be the sideways oscillatory 
motion of the vortices. During such sideways motion, dif- 
ferent portions of the vortices could cross the hot wire 
anemometer causing the fluctuations. These fluctuations 
were present during the measurements that led to Figure 28. 
This figure represents a uniform turbulence level across the 
channel. In Figures 29 and 30, the level of the turbulence 
was larger than in the previous cases. The magnitude of 
the fluctuations was also doubled or larger. The occasion- 
ally high value of the turbulence intensity would indicate 
turbulent flow. In other words, there was bursting to 
turbulent flow at the higher flow rates. 

The growth and motion of the vortices within the channel 
has been indicated by the mean velocity profiles and turbu- 
lence levels presented here. Evidence for the transition 
to turbulent flow has also appeared. 
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Figure 28. Turbulence Intensity Re = 1457. 
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Figure 29. Turbulence Intensity Re = 1722. 
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Figure 30. Turbulence Intensity Re = 1934. 



V. PRESSURE DROP MEASUREMENTS 



A. EXPERIMENTAL PROCEDURES 

To measure the pressure drop along the channel required 
a great deal of care. Each of the 28 pressure taps were 
cleaned and checked for tightness with high pressure air. 

The flow rate in the channel was adjusted and allowed to 
stabilize for approximately one hour before data was taken. 
The micromanometer was positioned on the table, leveled, 
and filled with the special fluid provided by the manufac- 
turer. 

The meniscus of this instrument was viewed under a 
magnifying glass and carefully adjusted so that the curve of 
the meniscus just touched the hair line on the glass tube. 
This was established as the zero position. It was found 
that the movement of even a small weight on the laboratory 
table would effect the level and therefore the zero of this 
micromanometer. Once the micromanometer was zeroed, no 
weight movement was permitted and the investigator did not 
touch the table while measurements were taken. 

The pressure fittings for each of the pressure taps were 
numbered and taped to the edge of the channel. The pressure 
drop between pairs of pressure taps at one foot intervals 
was measurable, but contained too much uncertainty to be of 
use. Typical values of 0.007 inches of water were measured 
as the pressure drop per foot. However, values of 0.001 
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inches higher or lower could be obtained in the same posi- 
tion of the channel at the same flow rate. It is possible 
that part of these variations were due to small fluctuations 
in the flow rate caused by the centrifical blower. Other 
causes were considered, but such fluctuations could not be 
positively identified or eliminated. 

The difference between the local pressure at each tap 
and atmospheric pressure proved to be more useful. At 
successive locations down the channel the magnitude of this 
pressure difference became larger and hence the relative 
error was less. In taking each pressure reading, the zero 
of the manometer was checked and then the flexible tubing 
from the micromanometer was connected to the pressure fit-' 
ting. The pressure tap number was checked and the inlet 
temperature reading was taken and both were recorded. The 
manometer reservoir level was adjusted with the fine screw 
dial until the meniscus of the fluid had returned to its 
zero location. The rotometer was checked to insure that 
the flow rate was correct and stable. The meniscus of the 
fluid was then rechecked and if it had remained in the zero 
position, the pressure indicated by the dial was recorded 
along with that flow rate. The micromanometer was discon- 
nected from the pressure fitting and then adjusted to 
return the meniscus to its zero position. If the dial 
indication was again zero, the pressure difference obtained 
was considered to be correct. Occasionally the micromano- 
meter would not re-zero and the pressure measurement at 
that location was retaken after the instrument was adjusted. 
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The pressure drop over two foot portions of both the 
straight section and the curved section were taken to verify 
the other results. For measurement over these intervals 
the upstream pressure tap was connected to the top of the 
micromanometer. The downstream tap was connected to the 
reservoir of the micromanometer as in the previous setup. 

B. PRESSURE DROP 

The pressures measured were plotted downward on Figure 
31. The top of this graph represents atmospheric pressure, 
and the amount that the static pressure in the channel was 
below the atmospheric level can be seen. The pressures are 
plotted against pressure tap number. The 28 pressure taps 
were three inches apart along the channel, and numbers one 
through 16 were on the straight portion of the channel. 

The pressure for six different Reynolds numbers were plotted, 
and the same patterns can be seen in each. The irregulari- 
ties in the curves around pressure taps 14 and 15 and around 
taps 18 and 19 were of a smaller magnitude than the uncer- 
tainty for these measurements. Pressure taps 14 and 15 were 
in the straight portion of the channel and no cause for the 
irregularities were apparent. However, burs or other dis- 
turbances which were undetected during construction of the 
channel may have caused these irregularities. Pressure taps 
18 and 19 were in the curved portion of the channel. In 
this case, the cause for the pressure irregularities may 
have been a difference in the position of the pressure taps 
relative to the vortices. A pressure tap that was being 
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Figure 31. Pressure Below Atmospheric vs. Pressure Taps. 



affected by the stagnation of outward moving fluid and one 
under the influence of moving fluid would probably give dif- 
ferent pressure indications. Further study would be required 
to confirm that this was the cause of the irregularities. 

The slope of the pressure curve in the straight section 
should compare with the exact analytic solution for laminar 
flow between parallel flat plates. The results of such an 
analysis [Ref. 6] show that the pressure drop is 

dp _ , yU 

dx 2 

a 

where a is half the plate spacing and y is the viscosity of 
the fluid. U is the mean velocity of the parabolic velocity 
profile between the parallel flat plates. However, the 
experimentally determined U was the volumetric flow rate 
through the rectangular channel divided by the cross sec- 
tional area. For a rectangular channel of aspect ratio 40, 
the difference in these two values of U would be small. 

The analytic value of U for a rectangular channel could be 
obtained from an infinite series given in the literature 
[Ref. 19]. The first term of this series would reduce U 
for the rectangular channel to 98.5315 percent of the value 
for an infinite parallel plate construction. The remaining 
terms of this series would rapidly become smaller. The 
experimental uncertainty was larger than this analytic 
difference, so that no correction was applied to the theo- 
retical pressure drop. 
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The experimental pressure drop was calculated by a com- 
puter program for both the straight and the curved portions. 
This program fit a least square straight line to nine con- 
secutive pressures representing a two foot section of the 
channel. The straight section of the channel was represented 
by pressure taps three through eleven inclusive. This 
avoided both the entrance region and the irregularities that 
were noted before. The curved portion was represented by 
pressure taps 19 through 27, where the vortices were most 
developed. The slope of each line was considered as the 
pressure drop in that portion of the channel. This program 
and its results are shown in Appendix C. These pressure 
drops in inches of water per foot of length, along with the 
analytic pressure drop calculated with the experimental 
values of U are shown in Table I. The percentage increase 
in the pressure drop for the curved section was calculated 
with respect to both the theoretical value and .the measured 
value of pressure in the straight section. These results 
were also tabulated in Table I. The negative values shown 
were probably due to experimental uncertainties. 

Because more significance might be obtained from a com- 
parison of pressure drops in the straight section with 
pressure drops in the curved section, where the vortices 
were fully developed, the same computer program was used to 
consider one foot section of the channel. The straight 
portion was represented by pressure taps six through ten, 
and the curved section by the last one foot, pressure taps 
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PRESSURE DROP PER FOOT FOR TWO FOOT INTERVALS CALCULATED BY LEAST SQUARES 



24 through 28. The computer results for these calculations 
are shown in Appendix C. The resulting pressure drops and 
percentage increases are shown in Table II. These results 
showed no negative values for the increase, and had the same 
magnitude as the other results. 

To verify and compare the other results, direct measure- 
ments of the pressure drop in the two foot intervals used 
above were made. The error in this type of measurement could 
have a magnitude as high as 0.002 inches of water. Never- 
theless, the results, along with the resulting percentage 
increases in the pressure drop, are presented in Table III. 
The values of the pressure drop measured in the straight 
section of the channel compared moderately well with the 
theoretical values. This comparison further supported the 
methods used to obtain the data. 

The percentage increase in the pressure drop from the 
straight section to the curved section was arrived at with 
three different measurements , and calculated with both the 
theoretical and the measured values in each case. The six 
different sets of percentage increases were plotted against 
the Reynolds numbers in Figure 32. The average increase for 
each Reynolds number is also identified in Figure 32. These 
results have considerable scatter due to the effects of 
experimental uncertainty. With an increase in the Reynolds 
number, there was an increased influence of the Taylor- 
Goertler vortices on the pressure drop. Even though the 
secondary flow was present, the increase in pressure drop in 
the curved section was small at the lower flow rates. At 
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PRESSURE DROP PER FOOT FOR ONE FOOT INTERVALS CALCULATED BY LEAST SQUARES 
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PRESSURE DROP PER FOOT FOR TWO FOOT INTERVALS MEASURED VALUES 
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higher Reynolds numbers, where the vortices became more 
active, the pressure drop in the curved section increased 
more than in the straight section. 

This data was nondimensionalized as the pressure drop 
coefficient, which was defined as 




where p was the pressure measured at each pressure tap. 

There was a larger error bound on these values because the 
uncertainty in U was added to that in the pressure drop. 

The calculation of the pressure drop coefficient was done by 
hand on a desk calculator. A sample calculation for the 
lowest Reynolds number is given in Appendix C. The results 
were plotted in Figure 33. 



83 




q.U0TOTg:g:aoD do^a aanssejcj 



cn 

a 

ftf 

Eh 



CD 

U 

2 

CO 

CO 

(U 

Jh 

a. 



84 , 



Figure 33. Pressure Drop Coefficient vs. Pressure Taps . 



VI . HEAT TRANSFER 



A. EXPERIMENTAL PROCEDURES 

In order to determine the effect of Taylor-Goertler 
vortices on heat transfer it was necessary to measure the 
heat transferred from a uniform temperature surface in one 
wall of the straight portion of the rectangular channel. 

The heat transfer coefficients of interest were defined by 
the equation 

Q = H * A* AT 

where Q was the heat convected, H was the heat transfer 
coefficient and A was the area of the heated surface. The 
temperature difference, AT, used depended on the particular 
heat transfer coefficient to be found. For the overall 
heat transfer coefficient, the log mean temperature differ- 
ence was used. The mean heat transfer coefficient was 
defined with the difference equal to the wall temperature 
minus the inlet bulk temperature of the fluid. The temper- 
ature difference used for the local heat transfer coefficient 
was the difference between the temperature of the individual 
heating element and the bulk temperature of the fluid enter- 
ing the control volume below that element. The differences 
in these definitions and the areas used will become clear 
during the discussion of the results. 

A heater plate with 16 separate resistance heating 
elements placed along the downstream length of the heater 
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was designed and built as described earlier. To measure the 
heat transferred from a uniform temperature surface it was 
necessary to control the voltage to each of the 16 heating 
elements. The resistance and voltage at each heating strip 
was to be measured to determine the power generated. The 
temperatures were to be measured at two locations on each 
heating elements. 

After the heater was in place in the channel wall, the 
resistance of each heating elements was measured with a 
Rosemont Commutating Bridge, Model 920A. The high current, 
regulated D.C. power supply fed the power controller which 
governed the voltage across each heating element. This 
specially constructed equipment had never been previously 
used or tested. The output level of each of the 20 circuits 
was set as low as possible. The input voltage was adjusted 
until the power controller's reference voltage was between 
the maximum recommended value of 8.800 volts and 8.700 volts. 
To insure that the heater plate did not become too hot, flow 
was established in the channel. The voltage in each of the 
first 16 power circuits was slowly increased, and the tem- 
perature along the plate was checked frequently. At any 
given power setting, four to five hours were required for 
the heated to come to equilibrium. As expected, in order 
to maintain the same temperature for each strip, the circuits 
near the leading edge of the plate required more power than 
the other circuits. If the variation of the power settings 
was not correct, the temperature distribution approached by 
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the plate was not uniform. When it was determined that the 
temperature would not be uniform, the power to the low 
temperature heating elements was increased and another 
settling period was allowed. Considerable time was required 
to obtain the first uniform temperature distribution. There- 
after, to obtain higher temperatures across the plate, each 
potentiometer's adjustment screw was turned equal amounts. 

As equilibrium was approached, small adjustments were made 
to remove temperature irregularities along the plate. In 
this manner, a uniform temperature could be achieved in a 
minimum time of about six hours. 

During continued testing of the equipment, it was found 
that the heater plate could sustain temperatures as high as 
200° F without damage. Higher operating temperatures were not 
attempted. During the early testing of the power controller 
the power circuit voltages were found to fluctuate even 
when the reference voltage was steady. Examination of the 
power circuit voltage showed a 60 cycle A.C. voltage super- 
imposed on the controlled D.C. voltage. The source of this 
A.C. voltage was the grounded side of the regulated power 
supply. Consequently, downstream from the protective 
diodes, a 0.05 yf capacitive connection was made across the 
inputs to the power controller. This connection removed the 
A.C. voltage and the fluctuation in the power levels. When 
the power circuits were loaded, the stability was found to 
be approximately ±.004 volts over a 30 minute period. The 
power available was actually more than required for these 
tests . 
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After the equipment was tested, data collecting began. 

At the beginning of each data run, the flow rate was estab- 
lished and power was supplied to the heater. After two 
hours, an ice bath was made and the temperatures along the 
plate were periodically checked. During the third hour, 
adjustments were made in the temperature and flow rate. The 
reference voltage was also checked at this time so that its 
variation could be determined later. 

The temperatures measured and recorded included the bulk 
temperatures of the entering and exiting air, the tempera- 
tures at the 32 thermocouple locations along the heater, and 
the temperatures of the four locations on the outside of 
the heater plate. The thermocouples were connected to a 
digital voltmeter by a set of switches. The switching and 
recording of the thermocouple outputs were accomplished by 
hand. As equilibrium was approached, the temperatures were 
recorded at 30 minute intervals. When the temperatures had 
stopped changing within the limits of accuracy of this 
experiment, the values were recorded as data. The reference 
voltage and flow rate were rechecked to insure no change had 
taken place. Then the digital voltmeter was connected to 
the power controller and the voltage across each heating 
element was recorded. The digital voltmeter was reconnected 
to the switches. The temperatures were recorded one addi- 
tional time, and the flow rate was again checked. The data 
was considered valid if no changes had taken place. The 
process of taking data took approximately 30 minutes. 
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B. RESULTS 

In order to calculate heat transfer coefficients from 
the data collected, it was necessary to account for the heat 
losses. In the geometry of this experiment, there were 
several potential losses. The largest was expected to be 
the conduction through the back side of the heater plate. 

For this reason, the insulation, described earlier was placed 
over the heated area. To obtain an estimate of the conduc- 
tion losses through the plate, a finite difference model of 
the experimental situation was used. 

A computer program was used to obtain a numerical solu- 
tion of the transient and steady state temperature distribu- 
tions and the heat flows of the three dimensional model. 

The program used is entitled TRUMP and was developed at the 
Lawrence Radiation Laboratory [Ref. 23]. This program was 
implemented for the IBM 360 computer at the Naval Postgrad- 
uate School by C. Erbayram [Ref. 24]. Because .of the length 
of this program and its output, only the steady state results 
of one run are included in Appendix D. 

In the computer program, one half of the symmetric 
arrangement of the heater plate, the insulation, and the 
surroundings were modeled by a total of 197 node points and 
606 thermal connections. The initial temperature of all 
the nodes, the heat generated at the nodes representing the 
heating elements, and the heat transfer coefficients repre- 
senting convection to the surroundings, were input data. 

The results were the steady state temperature at, and the 
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heat flow from each of the nodes. The time to reach steady 
state was also given by the TRUMP program, and compared well 
with the six hours required in the actual experiment. 

After the TRUMP program was tested, the heat generation 
rates actually measured were input to the program. The heat 
transfer coefficients representing the convection to the air 
in the channel were initially estimated and then iterated 
until the wall temperature arrived at by the computer program 
was the same as that actually measured. In this way, the 
TRUMP program was used to model each experimental run. The 
heat flows from the nodes representing the heating elements 
were used to calculate losses in the reduction of the exper- 
imental data. This process will be described next. 

Temperatures and several heat transfer coefficients were 
computed from the data collected by the computer program in 
Appendix D. A table of the thermocouple calibration data 
was included in the computer program. All of the tempera- 
tures were calculated by interpolating the thermocouple 
readings in this table. The local temperatures along the 
wall and the average wall temperatures are shown in the 
computer results contained in Appendix D. 

The log mean temperature difference was calculated from 
the equation 



TJim 



Tbo - Tbi 



Tw - Tbi 
Tw - Tbo 



where Tw was the average wall temperature, and Tbi and Tbo 
were the bulk temperatures of the air into and out of the 
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test section respectively. Using a table of fluid prop- 
erties vs. temperatures, the fluid properties were evaluated 
at the average bulk temperature. The overall heat transfer 
coefficient was then calculated from the equation 



the heated surface. The mean heat transfer coefficient was 
also calculated in the program. With the mean temperature 
difference defined as 



the mean heat transfer coefficient was calculated from the 
equation 



The symbols are the same as used for the overall heat trans- 
fer coefficient. These two values appear in the computer 
output. 

To calculate the local heat transfer coefficients, the 
power generated at each heating element was calculated from 
the resistance of the strip and the voltage across it. Two 
values of heat flow obtained from the TRUMP program were 
also read into the computer for each heating element. Since 
one node in the TRUMP program represented two adjacent heat- 
ing elements, the heat flow values from the TRUMP program 
were divided in half to represent the losses from the indi- 
vidual heating elements. At the end elements, the conduction 




cp * V * (Tbo - Tbi) 
Aw * T£m 



where V is the volumetric flow rate, and Aw is the area of 



Tm = Tw - Tbi 




cp * V (Tbo - Tbi) 
Aw * Tm 
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losses along the channel wall were attributed entirely to 
the first and last heating element. The two heat flows for 
each heating element represented the conduction loss and 
the heat convection into the channel. A loss factor was 
calculated by the computer program by dividing the conduc- 
tion loss by the sum of the heat flows. The power generated 
in each strip was multiplied by one minus the loss factor 
to give the useful heat at each element. 

The radiation loss was accounted for by the equation 

Qr = e * c * As * (Tw 4 - Tbi 4 ) 

where e was the emissivity of the nichrome taken as 0.7 and 
where As was the surface area of each heating element, 
0.014783 square feet. The temperatures used were converted 
to degrees Rankine, and the Stef an-Boltzmann Constant, o, 
was equal to 0.199967 x 10 _1 ° BTU/min ft 2 °R. The radiation 
loss was subtracted from the heat generated at each element 
less the heat conducted out of the plate to give the heat 
convected into the channel. 

The last loss accounted for was the convection from the 
heated air in the channel to the surroundings. Two thermal 
resistances were calculated for this type of loss. The 
first included terms for the forced convection in the chan- 
nel, the conduction through the plexiglas lower wall, and 
the free convection from the horizontal outside surface of 
the plexiglas. The second thermal resistance included terms 
for the forced convection in the channel, the conduction 
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through the aluminum side walls, and the free convection 
from the vertical outside surfaces. The forced convection 
was represented by the mean heat transfer coefficient 
calculated by the program. The free convection was repre- 
sented by correlation formulas given in the text, Heat 
Transfer , by J. Holman [Ref. 20]. The convection loss from 
the air through the walls of the channel was calculated from 

_ A * (Tb - 70.0) 
y c - 

where Tb was the bulk temperature of air entering the con- 
trol volume below each heating element, R was the thermal 
resistance for each case, and A was the corresponding area. 
The total loss of this type was the sum of the two Qc's 
calculated with each of the two thermal resistances. In 
summary the heat convected into the channel from each heat- 
ing element, designated as p(j), was the power generated 
at each element minus the conduction losses through the 
plate and the radiation losses. The heat which raised the 
bulk temperature of the air, was this convected heat, p(j), 
minus the convection losses from the air. 

The new bulk temperature of the fluid after passing a 
heating element was therefore calculated by the computer 
according to 

Tb ( j + 1 ) = Tb(j) + tP . 

J J p* cp* v 

The local heat transfer coefficient was then easily calcu- 
lated from the equation 
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Hf-n = p* cp* V * [Tb (j + 1) - Tb ( j ) ] 

J ' As * [Tw - Tb ( j ) ] 

where As was the area of one heating element. The distance 
of each heating element from the leading edge of the plate 
was also calculated and printed out along with the power 
convected, the bulk temperature of the fluid, and the local 
heat transfer coefficient. The average of these local heat 
transfer coefficients was also displayed on the computer 
outputs found in Appendix D. 

Five data runs were taken and reduced in the manner 
described above. The five runs cover four different Reynolds 
numbers and a moderate range of wall temperatures. The 
resulting local heat transfer coefficients were plotted 
against the heating element numbers in Figures 34 through 38. 
The center of the first heating element was 0.137 inches 
from the leading edge of the plate, and each following ele- 
ment was 0.274 inches farther downstream. 

The local heat transfer coefficients tend to be large 
at the beginning of the plate, decreasing more slowly with 
increasing distance along the heater plate. There is con- 
siderable scatter on these graphs due to the experimental 
uncertainty . 

The heat transfer coefficient at the second heating 
element appeared consistantly low. This was due primarily 
to the method of calculating loss factors, which did not 
charge this heating element for any of the end losses. How- 
ever, insufficient information about the heat flow was 
available to improve the technique. Other factors affecting 
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Figure 34 . Heat Transfer Coefficient Tw = 125 .- 




(Jo 7 ^3 u-rui/nia) guaTorggaoo aajsupa^ geeq 

c 



96 



Figure 35. Heat Transfer Coefficient Tw = 133. 
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Figure 36. Heat Transfer Coefficient Tw = 143. 
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Figure 37. Heat Transfer Coefficient Tw = 139. 
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Figure 38. Heat Transfer Coefficient Tw = 131. 



the scatter included the variation of temperature along the 
wall. This variation was as high as ±2.0° F and was probably 
due to small irregularities in the surface of the plate. 

The first nichrome strip on the fifth heating element was 
known to be indented, and its temperature was usually 1.5 
to 2.0 degrees higher than the surrounding nichrome ribbons. 

Typical error bounds are included in each graph. Dis- 
pite the scatter, the decrease in the local heat transfer 
coefficient is very much like that for a flat plate. How- 
ever, the heat transfer coefficients are consistently higher 
than those for a flat plate. 

A brief attempt was made to find some result in the 
literature with which to compare the heat transfer coef- 
ficients. The majority of the literature for developing 
heat transfer in rectangular channels considers both walls 
being heated. The Nusselt number of thermally fully 
developed flow between parallel flat plates with one wall 
at constant temperature was found to be 4.0 [Ref. 16]. For 
this experiment, a Nusselt number of 4.0, based on hydraulic 
diameter, corresponds to a heat transfer coefficient of 
0.0244 BTU/min ft 2 °F. For the lowest Reynolds number 
tested, the local heat transfer approaches this value. For 
the other Reynolds numbers, the local heat transfer coeffi- 
cient is higher than the fully developed value. 

This data was not nondimensionalized because it was not 
compared to any other results. The intended use of this 
data was for comparison with similar heated plate data from 
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the curved section of the channel. The procedures used to 
obtain these results were considered adequate for comparison 
to later work. However, several possible improvements were 
identified by this testing. Such improvements would include 
designing a system with fewer heating losses and developing 
a construction technique with would result in a smoother 
heated surface. 
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VII. CONCLUSIONS 



The observations made during this investigation have led 
to several conclusions about the nature of Tay lor-Goertler 
vortices in a curved rectangular channel. The vortices were 
observed to develop with distance along the channel. The 
importance of this observation was that an analytic descrip- 
tion of the flow would have to account for this growth. The 
vortices were observed to increase in amplitude and complex- 
ity with increasing flow rate. Thus at a higher Reynolds 
number, less distance along the channel was required to 
obtain the same degree of development. Several identifiable 
patterns of the developing flow were observed. The patterns 
were useful in determining the degree of development of the 
vortices . 

The study of the velocity profiles in the channel indi- 
cated that Taylor-Goertler vortices are a complex laminar 
flow. Furthermore, at flow rates below that necessary to 
sustain turbulent flow, momentary bursting to turbulent 
flow was observed. 

The effect of Taylor-Goertler vortices on the pressure 
drop in the curved rectangular channel was observed to vary 
from no effect to approximately a 20 percent increase. The 
amount of the increase in pressure drop due to the Taylor- 
Goertler vortices increased with increasing Reynolds numbers. 
This seems reasonable considering that as the secondary flow 
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increases in activity the interaction of the Reynolds 
stresses with the mean flow intensifies. 

The investigation of heat transfer from a constant 
temperature plate in one wall of the channel resulted in a 
procedure that would be adequate for determining the effects 
of Taylor-Goertler vortices on heat transfer. The heat 
transfer coefficients obtained for the straight section of 
the channel showed the trends expected and compared with 
the limiting data available. Improvements in the heat 
transfer measuring procedures could be made. 
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VIII. RECOMMENDATIONS 



At the end of this investigation, it was clear that 
several recommendations could be made regarding continued 
study of Taylor-Goertler vortices in a curved rectangular 
channel. More flow visualization should be done to clarify 
the patterns in the vortex flow. Different aerosol injec- 
tion techniques, including a wall slot, should be tried. 

More intense lighting and perhaps narrower slits should be 
used to obtain clearer pictures of the flow. An improved 
flow visualization technique should then be used to examine 
more locations along the curved portion of the channel. Due 
to the fact that there were motions in the flow patterns, 
still photography could not indicate all of the observations. 
A motion picture study of flow patterns would be ideal. 

More could be learned from a hot wire anemometer study 
where the probe could be moved sideways across the channel 
as well as radially in and out. In this manner, the perio- 
dic variation in the spanwise direction could be determined. 
The size of the vortices and their locations could then be 
determined. In addition, by knowing where the probe was 
in a vortex, the magnitude of the other velocity components 
could be determined. 

With respect to pressure drop, the placement of more 
pressure taps along the last portion of the straight section 
and along the entire curved section would improve the 
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reliability of the average pressure drop obtained. An 
investigation should be made to determine if the pressure 
variation in the direction perpendicular to the flow could 
be measured. 

It was clear that work should continue towards deter- 
mining the effect of Taylor-Goertler vortices on heat trans- 
fer. A heater plate for the curved section could be built 
in the same manner as the flat heater plate used in this 
experiment. However, it is recommended that a casting 
material with less shrinkage and better thermal insulation 
properties be used. If this is done, a smoother heated 
surface and fewer thermal losses would result. It is 
recommended that the development of a good heat transfer 
measuring technique be given considerable thought. 

For any analytic solution attempted, it was felt that 
it should account for the development of the vortices with 
distance and flow rate. If the magnitudes of the velocity 
components are determined, as suggested above, they should 
compare with any analytic solution obtained. The inter- 
action of the Reynolds stresses with the mean flow would 
have to be accounted for to predict the increase in pres- 
sure drop. Only a non-linear analysis such as J. Stuart's 
work [Ref. 25] could be expected to give such results. 

Other analytic techniques including numerical methods ought 
to be tried. 
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APPENDIX A 



ERROR ANALYSIS 



A calculation of the uncertainty for each of the major 
variables of this experiment was made in accordance with 
the methods of S. Kline and F. McClintoch [Ref. 26]. The 
estimates of the uncertainty in the measured quantities were 
made quite conservative so that there was considerable con- 
fidence in the uncertainties arrived at. As an example of 
the calculations the uncertainty in Reynolds number is 
calculated below. The Reynolds number is given by 



and the uncertainty is taken as 



dRe _ / dU 2 dDh 2 dv 2 

Re V U Dh v 

The uncertainty in' the flow velocity is obtained from other 
calculations like this one based on estimates of the uncer- 
tainty in reading the volumetric flow rate through the roto- 

meter and the uncertainty in the cross sectional area of the 

dU 

channel. The results of those calculations were -y- = 0.02805. 
The uncertainty in hydraulic diameter and kinematic visco- 
sity are obtained from estimates as 0.00426 and 0.00136 
respectively. Then the uncertainty in the Reynolds number 
is 
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dRe 

Re 



0.02840 



= J (7.868 + 0.1815 + 0.0185) x 10 4 = 



Other quantities and their uncertainties are given below 



The error bound for any quantity can be 


obtained by multi- 



plying its magnitude by its uncertainty. 



Quantity 


Uncertainty 


Ac 


0.00412 


Aw 


0.00320 


Cp 


0.07515 


cp 


0.00010 


De 


0.02838 


dH 


0.00426 


H 


0.27150 


Ho 


0.03517 


p (pressure) 


0.05000 


Ph (heat generated) 


0.00502 


p(j) (heat convected) 


0.25502 


Re 


0.02840 


T 


0.00750 


Tfm 


0.02000 


Tw 


0.01667 


U 


0.02805 


V 


0.02775 


p 


0.00080 
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APPENDIX B 



COMPUTER RESULTS FOR VELOCITY AND TURBULENCE 

The computer program for reducing the hot wire anemometer 
calibration data is shown below. This program was adopted 
from a program written by Mr. Loyd Smith while he was study- 
ing at the Naval Postgraduate School. The resulting cali- 
bration curve is also shown. 

The hot wire anemometer data reduction program was also 
adapted from a program by Mr. Smith and is shown below. The 
reduced data for the seven velocity and turbulence profiles 
is then included. 
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DIMENSION PO ( 20 ) f OATE ( 20 ) 

DIMENSION X( 16 ) , Y ( 16 ) 

DIMENSION XG ( 6 ) , VG 2 ( 6 ) 

REALMS LA6ELI/8H / 

REAL LAbEL2/4H / 

REALMS I T 1 TLE ( 1 2 ) 

DATA XG/O .0 , 1 .0 ,2 .0,3 .0,4.0 , 5.0/ 

DATA SUMX , SUMX2 , SUMY , SGMXY/ 0 .0 , C . C , 0 . 0 , 0 . C/ 

DATA A/3.55447/ 

C READS IN PROBE DATA AND CALI BRAT I CN DATA 

READ! 5, 10CC) PD 
R E AD ( 5 , 1000 ) DAT E 

1000 FORMAT ( 20 A4 ) 

WRITE ( 6 ,2000) PD 

200C FORMAT (■ 1' PROBE DATA: • ,20A4) 

WRI TE (6 ,2001 ) DATE 

2001 FORMAT ( ' 1 , 1 DAT E : *,20A4) 

P E AD ( 5 , 1001 )N 

10C1 FORMAT ( 110) 

WR ITE (6,2002)A,N 

2002 FORMAT ( ' • , ' FLOW CONSTANT: • , 8X , F 10 . 3 , / , ' NLMBER OF DA 

l T A POINTS: ' , 1 1 0 , / / , 6X , ' H ' , 5X , 1 V • , 9X , ‘ U • , 9 X , ' X ' , 9 X , 

2' V ,// ) 

DC 10 1=1, N 

C READS IN EXPERIMENTAL VALUES OF TP E MANOMETER READING 

C AND THE DC VOLTAGE OUTPUT OF THE ANEMOMETER. THE CCN- 

C STANT A IS SET FOR h IN IN. OF P2C ON A STANDARD CAY 

READ! 5, 10 C 2 ) H , V 

1002 FORMAT ( 2F10. 3 ) 

U-A*SQR T ( H ) 

X ( I )=SQRT ( U) 

Y ( 1 ) = V**2 
X2=X( I ) **2 
X Y=X ( I ) *Y ( I ) 

SLMX=SUMX+X( I ) 

S L M X2 = S CM X 2+ X 2 
SCMY=SUMY + Y ( I ) 

SLMXY=SUMXY+XY 
WRITE(6,20C3)H, V,U,X( I ) ,Y(I) 

2003 FORMAT ( 5F 1 0 .3 , / ) 

1C CONTINUE 

C COMPUTE CONSTANTS BY LEAST SQUARES ANALYSIS 



D=N*SUMX2- ( SUMX- *2 ) 

BETA=( N*SUMXY-SUMX*SUMY) /D 
VC2=( SUMY *SUMX2- SUMXY -SUMX ) /D 

WRITE(b,2004)6ETA,VC2 , 

2004 FORMAT ( ' ','bETA = ' , F 1 0 . 3 , 1 OX , ' VC2 = • , F 1 C . 3 ) 

RE AD (5, 1003) ( ITITLEI I) ,1 = 1,6) 

READ! 5, 10C3 ) ( IT I TL E ( I ) , 1 = 7 , 12 ) 

1003 FORMAT ( 6A 3 ) 

C POINT PLOTS DATA , , „ „ , 

CALL DR A W ( N , X , Y , 1 ,2 ,LA6EL1 , ITITLE, 1.0, 1 .0 , 1, 1,0, 0,4, 

C 3 fPE’LEAsl SQUARES EQUATION IS USEC TO COMPUTE SIX 
C POINTS FOR A STRIGHT LINE 



DC 20 1=1,6 , , 

VG2 ( I ) = V02 + BETA-XG ( I ) 



^ CALL*DRAW(o,XG,VG2,3,0,LABEL2, ITITLE, 1.0, l.C, 1, 1,0,0, 
44 ,5,0, LAST) 

S TCP 
END 
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Voltage squared 




Hot wire anemometer calibration curve 
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DIMENSION DAT A L( 20) , CATA2 ( 20 ) , H ( 25 ) , U ( 25 ) ,7(25) 

REAL* 8 ITITL1 (12) , ITITL2(12) 

REAL LABEL/4H / 

30 READ(5 , SOCJNPROB 
SCO FORMAT ( I 1C) 

I F ( RlPRO B . L T . 0 ) STOP 
REAC(5, lOCCJDATAl 
R E AD ( 5 , 1000 ) DAT A2 
1000 FORMAT ( 20 A4 ) 

WRITE (6 ,2000) 

2000 FORMAT (• 1 ',////////// ) 

WRITE(6 ,2001 ) D AT A 1 
WRITEIo ,2001 JDATA2 
20C 1 FORMAT ( ' ',21X,20A4) 

REAC(5, 1001) RO , VC2 , BET A 
1CC1 FORMAT ( 3F 1C . 6 ) 

WR ITE ( e> , 2002 ) RO,VC2,BETA 

2002 FORMAT (33X,' IMTICAL PCS I T I CN= • , F7 . 3 , / , 24X , ' V02 = ', 
2F 10 .6, 9X, 'dETA = ' , F 10.6,// ,30X, * H 1 ,19X, ‘O' , 1 9 X , ' T ' ,//) 



1C 

1002 



R C F = 0 

P EAD ( 5 , 10C2JR, V,E 
FORMAT ( 2F10.6) 

IF (R.LT .0.0) GO TC 20 
R C P = NDP + 1 
I = RiCP 
Fi(I )=R-RO 

U ( I ) = ( ( (V**2)-VC2 )/EETA)**2 
T(I)=(4.C*V*E)/((V**2)-V02) 

W R 1 TE ( 6 , 2003 ) F ( I ) , U ( I ),T( I ) 

2003 FORMAT ( 21X,F10.6, 1CX, F10.6, 10X, F1C.6,/ ) 

GO TO 10 
WR ITE (6 ,2004 ) 

FORMAT ( //////////) 

R E AC ( 5 , 10 C3 ) ( ITITLl(I), 1=1,6) 

READ! 5, 1003) ( 1TITL1 ( I ) , 1 = 7, 12) 

! FORMAT ( 6A8 ) 

WR ITE (6 ,2005 ) 
i FORMAT! ///////////////) 

CALL DRAW (R|DP,U,H, 0,0, LABEL, ITITLl, 1.0,0.05,0,0,0,0, 
88 ,5 ,0 , LAST ) 

R EAD ( 5 , 1 C C 4 ) ( IT I T L 2 ( I ) ,1 = 1,6) 

R E AC ( 5 , 10 04) ( ITITL2( I ) , 1 = 7, 12) 

FORMAT ( 6A8 ) 

CALL DRAW(RiDP,T ,H,0,0,LABEL, ITITL2, 0.010 , .05,0,0,0,0, 
68, 5,0, LAST) 

GO TO 30 
ERD 



2C 

2CC4 



20C ! 



1CC4 
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CGNCAVE TO CONVEX WALL DISTANCE, VELCCITY AND TLRBULENCE 
FRC6E NC. TSI1279 DATE 5 FEB 73 23? FLOW RATE RE= 6C9 



V02 = 


INITICAL PCSI T ICN= 0. 
2 .020000 BETA 


27 C 

= C. 710000 


H 


U 


T 


C .C1CCOO 


0.613295 


0.031175 


C .C2C000 


0.620402 


0.032164 


C . C3CCC0 


0.5516C6 


0.028536 


C.C4C000 


1.536921 


0.023217 


C .C5CCCC 


1.958789 


0.020265 


C.C6CCC0 


2.504096 


0.017980 


C .07 COOO 


2.665120 


0.017536 


C .C8CC00 


2.90C582 


0.016941 


C .C8C000 


3.111880 


0. 016465 


C.1CCCC0 


3 .239269 


0.0159 17 


C .1 ICCCO 


3.332104 


0.015737 


C. I2C000 


3.407528 


0.015596 


C . I3CCC0 


3.426528 


0.015561 


C.14CCC0 


3 .407528 


0.015596 


C.I5CC00 


3.3321C4 


0.015737 


C . 1 6 CC CO 


3.202550 


0.015990 


C . I7C000 


3.C40466 


0.016621 


C . 18C0C0 


2.815131 


0.017149 


C .190000 


2.504096 


0.017980 
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CONCAVE TO CONVEX WALL CISTANCE, VELOCITY AND TURBULENCE 
FRCEE NC. T S I 1279 DATE 5 FEB 73 33? FLOW RATE RE= 875 

I N I T I CAL P C S I T IGN = C.270 

VG2 = 2.C20000 BETA = C. 7IOOOO 



H 


u 


T 


C .CICOOO 


0.701554 


0.03C957 


C . C2CCC0 


1 .258702 


0. 027440 


C.C3CCC0 


I .545134 


0.023128 


C .C4CC00 


2.665120 


0. C1565C 


C .C5CC0C 


3.111880 


0.017332 


c.cecoco 


3.580913 


0.016380 


c .C7CC00 


3.543310 


0.015227 


c.ceccco 


4 .153802 


0.014526 


c.cscooo 


4.370939 


0.014625 


C . ICCOCO 


4.545487 


0.014152 


C.liCOOO 


4.640411 


0; 0I4C43 


C.I2C000 


4 .686278 


0.013589 


C . I3CCOO 


4 .663301 


0. 014C 16 


0 . 14C000 


4.572118 


0.013877 


C . I5C00C 


4.455676 


0.0145 14 


0.16C000 


4.283268 


0.014746 


C. I7CC00 


4 .026724 


0.015111 


c.ieccco 


3.715356 


0.015865 


C . 1 9CC00 


3 .369696 


0. 0165C5 
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CGNCAVE TO CONVEX WALL DISTANCE, VELOCITY AND TURBULENCE 
FRCBE NO. T S I 127 9 DATE 5 FEB 72 4 4? FLOW RATE RE=1166 

INITICAL PC S I T I CN= 0.270 
VC 2 = 2.C20000 BETA = C. 710000 



H 


U 


T 


0 .01C000 


0.942583 


0.032475 


C.C2CCC0 


1.958789 


0.025157 


C . C3CC00 


2 .815131 


0.021C60 


C .C4C00C 


3.659631 


0.018945 


C .C5CCOO 


4 .153802 


0 .0175C0 


C .06CG00 


4.594827 


0.0 166 18 


C .C7CCC0 


5.015061 


0.015774 


C .C8C000 


5 .258296 


0.015C16 


C .C9CC00 


5.483244 


0. 014213 


C . 1CCCC0 


5.688062 


0.014111 


C . 1 10000 


5 .844732 


0.012514 


C . 12CCC0 


5.977282 


0.013398 


C . 13CC00 


6.003999 


0.012152 


C . 14CCCC 


5.95C615 


0.012421 


C.15CC00 


5.844732 


0.012514 


C.16C000 


5.636421 


0. 0 127C4 


C. 17CCC0 


5.282988 


0.014266 


C .180000 


4 .919745 


0.014929 


C .19CCC0 


4.459676 


0. 0155 15 
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CONCAVE TO CONVEX WALL DISTANCE, VELOCITY AND TUREULENCE 
FRCBE NC. T S I 1279 DATE 7 FEB 72 ASS FLOW RATE RE=12S8 
IMTICAL PCS IT ICN= 0.270 
V02 = 2.C20000 BETA = 0.7LOOOO 



H 


u 


T 


0 .OIOGOO 


l .320405 


0. 02C549 


C.C2CCC0 


2 . C 84258 


0.027238 


C.C3C000 


3 .600484 


0.024517 


C .C4C00C 


4.825512 


0. 02 1825 


C .C5CCOO 


5.407611 


0.02C888 


C .C6C000 


5.87 1078 


0. C17S67 


C.C7CCC0 


6.C577C6 


0.017772 


C .C8CCOO 


6 .111675 


0.017265 


C .CSCCCO 


6 . 2205C4 


0. 0167 12 


C. 1CCCC0 


6 .3C2920 


0.016151 


C . I 10000 


6.235582 


0. 015229 


C . 12CC00 


6.4S7612 


0.015128 


C.13CC00 


6.582447 


0.014623 


C. 14CC0C 


6 .667787 


0.014123 


C.15CC0C 


6.639285 


0.014156 


C. 16C000 


6.554168 


0 . 0 1 42 26 


C . 17CC00 


6.2C2920 


0.014440 


C .18CC00 


5.977282 


0.015184 


C .ISCCCO 


5.457553 


0.015726 
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CCNCAVE TO CONVEX WALL CISTANCE, VELOCITY ANC TIRBULENCE 
PRCBE NC. T SI 1275 DATE 6 FEB 73 55? FLOW R £7 E RE=1457 



VC2 = 


IN I T I CAL PC S I T IGN= C.27C 
2.C2C000 BETA = C 


.710000 


H 


U 


T 


C.CICCOO 


2.348951 


0.023578 


C . C2CCOO 


4.872481 


0.015326 


C .C3CG00 


6 . 138762 


0.017661 


C . C4CCC0 


6.811521 


0.014869 


C .05C00C 


7 .045468 


0.014679 


C .C6CCOO 


6.927876 


0. 014772 


C.C7CCOO 


6.869542 


0.014821 


C .C8COOO 


6.725045 


0.014541 


C . C9CCCO 


6.457812 


0.015128 


C.1CCCOO 


6 .441764 


0.015188 


C.I1CCCO 


6 .330523 


0.015290 


C . 12CCC0 


6 .3 C2920 


0.015216 


C . 120000 


6 .275269 


0. C145C3 


C.14CCCC 


6.247511 


0.014528 


C.15C000 


6 .247911 


0.014528 


C.16CCC0 


6 . 193172 


0.014528 


C .17CCC0 


6.020813 


0.014687 


C . 18C000 


5.766069 


0.014943 


C . 15CCC0 


5.357558 


0.0163C5 
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CCNCAVE TO CONVEX WALL CISTANCE, VELOCITY ANC TURBULENCE 
PROBE NG. T S I 1279 OATE 6 FEB 73 65? FLOW PATE RE=1722 

IN I T I CAL POSIT ICN= 0.270 
V02 = 2.C20000 BETA = 0.710000 



H 


U 


T 


C. 010000 


2.488311 


0. 05C629 


C.C2CCC0 


5. C87305 


0.047533 


C.C3C000 


6.358225 


0.0436 13 


C .C4CC0G 


7 .254210 


0.04 1479 


C .C5CCC0 


7.497396 


0.032772 


0 .C6CC00 


7.466739 


0.032823 


C.C7CCC0 


7.405624 


0.024694 


C .C8C000 


7.284319 


0 . 0207 07 


C .C9CCC0 


7 . 164276 


0.020838 


C . 1CCCC0 


7 . 104720 


0.02CSC4 


C . 1 ICOOO 


7 .015953 


0. C189C3 


C.12CCC0 


6.957165 


0.018S64 


C .130000 


6.927876 


0.016864 


C . 14CCCC 


6 . £98683 


0.016511 


C . 15CCCC 


6.811521 


0.016953 


C . 16C000 


6.696387 


0. 017 103 


C. 17CCC0 


6.554168 


0.019359 


C .180000 


6.330523 


0.024C27 


C . 15CCC0 


5.924045 


0.024647 
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CCNCAVE TC CONVEX WALL CISTANCE, VELOCITY ANC TLRBULENCE 
PRCBE NC. TSII279 DATE 6 FEB 73 72? FLOW RATE RE=1934 



VC2 = 


INIT ICAL 
2 . C20000 


POSITION^ C.27C 
BETA = C 


.710000 


H 




U 


T 


C.C1COOO 




2.425794 


0.095928 


C.C2C0C0 




4.778824 


0.097417 


C .020000 




6 .0208 13 


0 . 0 8 9 C 1 1 


C . C4 CC CO 




7 . 164276 


0.075C16 


C .05CC0C 




7 .466739 


0.061543 


C .C6CCC0 




7 . 62C902 


0.052928 


C.C7CCC0 




7 .589893 


0.052CC9 


C .C8C000 




7.528167 


0. C49C62 


C . C9CCC0 




7.284319 


0.041414 


C.1CC000 




7.254210 


0.04 1479 


C . 1 ICCCO 




7.405624 


0.037C41 


C.12CCC0 




7 .405624 


0.0288C9 


0 .120000 




7.466739 


0. C2872C 


C.14CCCC 




7.466739 


0.024617 


C . 1 5CCCC 




7.558975 


0.02C420 


C . 16C0C0 




7.558975 


0. C2C420 


C .17CCC0 




7 .466739 


0.02C514 


C . 18C000 




7.214522 


0.C2C675 


C. 19CCC0 




7 . C 1 59 53 


0. 0252C4 
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APPENDIX C 



CALCULATIONS AND COMPUTER RESULTS FOR PRESSURE DROP 

The computer program written to find the slope of the 
pressure curve is shown below. The results for six dif- 
ferent Reynolds numbers are shown. The first set of results 
represent nine pressure readings or two foot intervals. The 
second complete set represents selected one foot intervals. 
The last item shown in this appendix is a tabular calcula- 
tion of the pressure drop coefficient for the lowest Rey- 
nolds number tested. 
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OOOOOOOO 



CALCULATION OF THE LEAST SQUARES SLOPE FOR PRESSURE 
DRCP DATA. DATA READ IN IS NUMBER uF DATA POINTS, 
REYNOLDS NUMBER, POSITION OF PRESSURE TAP IN INCHES 
ANC PRESSURE IN INCHES OF WATER. 



DIMENSION TYPE! 2 ) 

WRITE (6 ,11 ) 

I 1 FORMAT ( ' 1« ,////) 

L = 0 

2 CO R E AD ( 5 , 10 ) N , RE , TY F E (1) , TY P E (2 } 

1C FORMAT ( 16, F6. 1.2A4) 

I F (N.LT .1 ) GO TO 3C0 

S L M X= 0 . C 

SUM Y= 0 . C 

SUMX2 = 0 .0 

SLMXY-0 . 0 

L = L + 1 

1FIL.GT .3 ) GO TO 310 
GO TO 210 

3 1C WR ITE ( 6, 1 1 I 

L = 1 

2 1 C WRITE (6, 20) RE 

20 FORMAT ( 1 ',31X, 'SLOPE OF PRESSURE DROP CURVE FCR *, 

1 ' REYNOLDS NUMBER = ' »F7.1,///,35X, 'POSITION (FT)',20X, 
2 • PRESSURE ( I N H2G ) ' ,//) 

C = N 

DC 100 1=1, N 
READ<5,2 2) T , Y 

22 FORMAT (2F10 .5 ) 

X = T/12 . C 

WR ITE ( 6 , 23 ) X , Y 

23 FORMAT (39X, F7 .5 ,25X,F7 .5) 

X2=X**2 

X > = X*Y 
SLMX=SUMX+X 
SLMY=SUM Y+ Y 
SLMX2= S UMX2+X2 
SUMXY=SUMXY+XY 
ICC CONTINUE 

C= (SUMX }**2~C*SUMX2 
A C= ( S UMX* SUMXY-S UMY *SUMX2 )/ C 
A 1= ( SUM X^ SUM Y-C * S L MXY ) /D 
W R I T E ( 6 , 24 ) TYPE ( 1 ) , TYPE ( 2 ) , A1 

24 FORMAT! bX ,// ,33X, • THE PRESSURE DROP IN THE ',2A4,IX, 

3 ' SECT ION IS' , F 12 . 6 , // / / ) 

GC TO 200 
3CC STOP 
ENC 
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SLCPE OF PRESSURE DROP CURVE FOR REYNCLCS NUMBER 



609 .0 



POSITION (FT) 



0.0 

0. 250CC 
0.500CC 
0.75000 
I. OOOCO 
1.25000 
1.50000 
1.75000 
2.00000 



ThE PRESSURE DROP IN 



SLCPE OF PRESSURE CROP 



POSITION (FT) 



PRESSURE (IN 



0 .01 ICO 
0.01200 
C. 01300 
0 .014CC 
0.01450 
C. C15C0 
0 .01650 
C. 0180C 
C.015C0 



STRIGHT SECTION IS 



RVE FOR REYNCLCS NUMBER 



PRESSURE (IN 



0.0 

0.25000 
0. 500CC 
0.75000 
1. OOOCO 
1.25000 
1.50000 
1. 75C0C 
2. OOOCO 



C. 03000 
0.03200 
0.03400 
C . C35CC 
0.03650 
C. 038C0 
C.C350C 
0.04000 
0.04 150 



THE PRESSURE DROP IN THE CURVED SECTION IS 



SLCPE OF PRESSURE 
POSITION (FT)' 



DROP CURVE FOR REYNOLDS NUMBER 

PRESSURE (IN 



0.0 

0.25000 
0. 500CC 
0. 750CC 
1.00000 
1. 25000 
1. 5000C 
1 .750C0 
2. OOCCC 



C. 014C0 
C .C16CC 
0 .01800 
0. C1550 
C.021GC 
0.02300 
C. C250C 
0 .027CO 
0.02900 



THE PRESSURE DRQP IN THE 



STRIGHT SECTION IS 



F 20 ) 

0.003867 

= 6GS.0 

H20 ) 

0 .005533 

= 875.0 

H2C ) 

0 .007367 



12 .1 



SLCPE OF PRESSURE DROP 



POSITION (FT) 



0.0 

0.25000 
0. 500CC 
0. 750C0 
I .00000 
1. 250CC 
1. 500CC 
1.75000 
2. OOOCC 



THE PRESSURE DROP IN THE 



CURVE FOR R EVNGLCS NUMBER 



PRESSURE (IN 



C. 042CC 
0 .04400 
0.04550 
C.047CC 
0 .045CC 
0.05000 
C. C52CC 
0.C540C 
0.05600 



CURVED SECTION IS 



SLOPE OF PRESSURE DROP CURVE FOR REYNOLDS NUMBER 



POSITION (FT) 



0.0 

0. 250C0 
0.50000 

0. 75CCC 

1. 00000 
1.25000 

1. 5000 C 
I. 750CC 

2.00000 



THE PRESSURE DROP IN THE 



PRESSURE (IN 



0.02400 
0. C260C 
C.028CC 
0 .03050 
0. 033C0 
0 .03500 
0.03700 
C. C3900 
0.04200 



STRIGHT SECTION IS 



SLCPE OF PRESSURE DROP 



POSITION (FT) 



0.0 

C. 25000 
0. 500CC 
0.75000 
1 . 00000 
1.25000 
1.50000 
1.75000 
2.000C0 



THE PRESSURE DROP IN 



RVE FOR REYNCLCS NUMBER 



PRESSURE (IN 



C.C59C0 
0.06200 
C. 06300 
C.066CC 
0 .06900 
C. 07 ICC 
0 .073C0 
0.07600 
C. 079C0 



CURVED SECTION IS 



= 875. C 

H20 ) 



0 .0068C0 

= 1166.0 
F2G ) 



0. CC89C0 

= 1166. C 
h20 ) 



0 .009800 



122 



SLCPE OF PRESSURE ORCP CURVE FOR REVNCLCS NUMBER 



1484 .0 



PCS I T I ON (FT) 



0.0 

0.25CC0 
0. 50GCC 
0 .75000 
1 . 00000 
1.25000 
1.50000 
1. 75C00 
2.000C0 



THE PRESSURE DROP IN THE 



PRESSURE (IN 



0.03200 
0.03500 
C. 03800 
C.C4CCC 
0.04300 
0. C46CC 
C . C49CG 
0 .05200 
0.0550C 



STRIGHT SECTION IS 



SLOPE OF PRESSURE DROP CURVE FOR REVNCLCS NUMBER 



POSITION (FT) 



PRESSURE (IN 



0.0 

0 .25000 
0.50000 
0. 75CCC 
1.00000 
1.25000 
1. 500CC 
1.75000 
2 . 00000 



0.0780C 
C. C82CC 
0 .08500 
C.C88CC 
0. 093C0 
0 .09500 
C.C980C 
C. 102CC 
0 . 10 2C0 



THE PRESSURE CROP IN THE 



CURVED SECTION IS 



SLCPE OF 



POSIT 



0 

0 

0 

0 

1 

1 

1 

1 

2 



PRESSURE DROP CURVE 



ION (FT) 



.0 

. 250CC 
. 5 OOC C 
.75000 
. OOOCC 
. 25000 
.50000 
. 750CC 
. OOOCO 



FOR REVNCLCS NUMBER 



PRESSURE (IN 



0 .03500 
0.04300 
C. 047CC 
0.050CC 
0.05300 
0. 057CC 
0 . 06 ICC 
0.06400 
C. 067CC 



THE PRESSURE DROP IN THE 



STRIGHT SECTION IS 



H20 ) 

0.011400 

= 1484.0 
H20 ) 

0.0134CC 

= 1748 . C 
F20 ) 

0.014000 
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SLCPE OF PRESSURE DROP CURVE FOR REYNOLDS NUPEER 



1748 .0 



POSITION (FT) 



0.0 

0. 250CC 
0. 5000C 
0.75000 
1. OOOCC 
1.25000 
1.50000 
1.75CC0 
2.00000 



THE PRESSURE CROP IN THE 



PRESSURE (IN F 20 ) 



0.C93CC 
0.09800 
C. 10CC0 
0 . 102CC 
C. 1 1000 
C. 112C0 
0 .11600 
C. 12100 
0. 1260C 



CURVED SECTION IS 0.016133 



SLCPE OF PRESSURE CROP CURVE 



POSITION (FT) 



0.0 

0.25000 
0. 50CCC 
0. 75000 
1 . 00000 
1. 25CCC 
l. 500CC 
1 . 75 OCC 
2. OOOCC 



FOR REYNOLDS NUMBER = 1987. C 



PRESSURE ( IN h20 ) 



C. 048CO 
0.C52CC 
0.05600 
C. C6CCC 
0.064CC 
0.06800 
C. C73CC 
C.076CC 
0.08000 



THE PRESSURE DROP IN THE 



STRIGHT SECTION IS 



0.016133 



SLOPE OF PRESSURE DROP CURVE 



POSITION (FT)' 



0 . 0 

0.2500C 
0.50000 
0. 75CCC 
1.00000 
1.25000 
1. 500CC 
1.75000 
2.00000 



FOR REYNCLCS NUPEER = 1987.0 



PRESSURE (IN P2C ) 



0.10800 
C. 1 14CC 
C . U7CC 
0.12100 
0. 129CC 
0 . 132CC 
0.13700 
0. 142C0 
0.147CC 



THE PRESSURE DROP IN THE 



CURVED SECTION IS 



C. C194CC 



SLOPE OF PRESSURE DROP CURVE FOR REVNCLCS NUMBER 



609 .0 



POSITION (FT) 



0 . 0 

0.25000 
0.50000 
0. 75CCC 
I. 00000 



THE PRESSURE DROP 



SLCPE OF PRESSURE 



POSITION (FT) 



0.0 

0.25000 
C. 50CCC 
0.75000 
1 .00000 



THE PRESSURE DROP 



SLCPE OF PRESSURE 



PCSITIGN (FT) 



0 .0 

0. 250C0 
0.5000C 
0.75000 

1. ooocc 



THE PRESSURE DROP 



SLOPE OF PRESSURE 



POSITION (FT) 



0.0 

0.25000 
0.50000 
0. 7500C 
L. 00000 



THE PRESSURE DROP 



PRESSURE (IN F2C) 



0.01400 
C. 0L450 
C .C150C 
0.01650 
0. 018C0 



IN THE STRIGHT SECTION IS 0.004000 



DROP CURVE FOR REVNCLCS NUMBER = 609.0 



PRESSURE (IN H2C ) 



C.C380C 
C.C39CC 
C .04000 
C. 04150 
0 .G42CC 



IN THE CURVED SECTION IS 0.0C5CCC 



DROP CURVE FOR REVNCLCS NUMBER = 875.0 



PRESSURE (IN H20 ) 



C. 01950 
0.02100 
C.C23CC 
0 .025CC 
0.02700 



IN THE STRIGHT SECTION IS 0.0076C0 



DROP CURVE FOR REVNCLCS NUMBER = 875.0 



PRESSURE (IN F2C ) 



0.05000 
C.C52CC 
C.C54CC 
0 .05600 
C. 056C0 



IN THE CURVED SECTICN IS 0.008000 
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1166 .0 



SLCPE OF PRESSURE DROP CURVE FOR REYNOLDS NUMBER 



PCS I T I ON (FT) 



0.0 

0.25000 
0.500CC 
0.75000 
1. OOOCC 



PRESSURE (IN 



C. 03050 
0.03300 
0.C3500 
C .037CC 
0 .03900 



THE PRESSURE DROP IN THE STRIGHT SECTION IS 



SLOPE OF PRESSURE DROP CURVE FOR REYNCLCS NUN E ER 



POSITION (FT) 



0.0 

0.25000 
0.50000 
0.75000 
l. OOOCO 



PRESSURE (IN 



0.07100 
C.073C0 
C.C76CC 
0 .07900 
0. C82C0 



THE PRESSURE DROP IN THE CURVED SECTICN IS 



SLCPE OF PRESSURE DROP CURVE FOR REYNCLCS NUMEER 



POSITION (FT) 



0 . 0 

0.25000 
0.50000 
0. 750CC 
1.00000 



THE PRESSURE DROP IN THE 



PRESSURE (IN 



0.04000 
C.043C-G 
0 .0460 C 
0.04SC0 
C. C52CC 



STRIGHT SECTICN IS 



SLCPE CF PRESSURE CROP CURVE FOR REYNCLCS NUMBER 



POSITION (FT) 



0.0 

0.25000 
C. 500CC 
0.75000 
1 . 00000 



THE PRESSURE CROP IN THE 



PRESSURE (IN 



C.095CC 
C .098CC 
0.10200 
C. 105CC 
0 . 10900 



CURVED SECTION IS 



H20 ) 



C.CC84CC 

= 1166.0 
F20 ) 



0.01 1200 

= 1484.0 
H2C ) 



C.0120C0 

= 1484.0 
h20 ) 



0.014CCO 
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SLCPE OF PRESSURE CROP CURVE FOR REVNCLCS NLMEER 



1748 .C 



POSITION (FT) 



0.0 

0.250GC 
0.500CC 
0. 75CCC 
1.00000 



THE PRESSURE DROP IN THE 



FFESSURE (IN 



C. 05000 
C. 052CC 
C .057C0 
0.06100 
C . C64CC 



STRIGHT SECTION IS 



SLCPE OF PRESSURE DROP CURVE FOR REVNCLCS NUMBER 



POSITION (FT) 



0.0 

0.25000 
0. 500 C C 
C. 75000 
1 .00000 



THE PRESSURE CROP IN THE 



PRESSURE (IN 



C. 112CC 
0 . 1I6CC 
0.12 100 
C. 126CC 
0 . 13C0C 



CURVED SECTION IS 



SLCPE OF PRESSURE DRCP CURVE FOR REYNCLCS NUMBER 



POSITION (FT) 



0 .0 

0. 250CC 
0.5000C 
0.75000 
1. OOCCC 



THE PRESSURE DRCP IN THE 



PRESSURE (IN 



C .06CCC 
C .06400 
C.068CC 
0 .072CC 
0.07600 



STRIGHT SECT ICN IS 



SLOPE OF PRESSURE DROP CURVE FOR REYNCLCS NUMBER 



PGSITION (FT) 



0.0 

0.25000 
0.50000 
0. 75000 
1.00000 



THE PRESSURE DROP IN THE 



FFESSURE (IN 



0.13200 
C. 137CC 
C.142C0 
0 .14700 
C. 15200 



CURVED SECTION IS 



F2C ) 

0.014400 

= 1748.0 
H20 ) 

C.C184CC 

= 1987 .0 
F20 ) 

C .C164CC 

= 1587 .0 
F 20 ) 

0 .020000 
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APPENDIX D 



COMPUTER RESULTS FOR HEAT TRANSFER 

The final results for one run of the TRUMP program are 
shown on pages 130 through 144. The computer program for 
the calculation of the heat transfer coefficients follows 
the TRUMP program. The five sets of heating data results 
are presented with the temperature distribution shown on 
the first page and the heat transfer coefficients on the 
second page of each of the results. 
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CALCULATION OF OVERALL AND LOCAL FEAT TRANSFER COEF- 
FICIENTS FROM DATA FOR A WALL HEATER IN A RECTANGULAR 
CHANNEL . 

THE DATA REQUIRED INCLUDES ThE HEATING STRIP RESIS- 
TANCE AND VOLTAoE DRGPS TO GIVE THE POWER, THE THERMO- 
COUPLE MILLIVOLT REACINGS AND CALCULATION DATA, ANC 
CONSTANT VALUES FOR AREA, FLOW RATE, FLUIC PROPERTIES 
TABLES, AND LOSS FACTOR CALCULATIONS. 

THE UNITS ARE FEET, MINUTE, BTU , LBM , AND CEGREES F. 



DIMENSION TRC( 14, A) ,R( 181 ,V( 18) , TM ( 40 ) , TD ( A C ) » TB ( 2 C ) 
DIMENSION h!18),U(13)»X(18),P(18),GT(18) 

DC 100 1=1,14 

READ (5, 10 I TRC( I, 1) ,TRC( I,2),TRC(I,3),TRC!I,4) 

1C FORMAT ! F10. 1,2F10.5,F10.3) 

ICC CONTINUE 

CC 200 J= 1 , 1 8 
P EAC ( 5 , 12 1 R ( J I , V ( J ) 

12 FORMAT! 2F 10.4) 

2CC CONTINUE 

DO 340 K = 1 ,40 
R EAC( 5 , 14) TM ( K ) 

34 C T C ( K ) =0 .0 
14 FORMAT! F10. 3) 

DC 320 L=.L , 20 
32C T E ( L ) =0 . 0 

DC 318 M= 1 , 18 
H (M ) = 0 .0 
318 P ( M ) =0 . C 

DO 322 L=l,18 
322 R EAC! 5, 20) U(L) , G T ( L ) 

2C FORMAT (2F8.3) 

READ! 5, 11) XO , DX » A N , AS 
11 FORMAT (4F10. 6) 

R E AC ( 3 , 13 ) u » TO , S , A 1 , A2 

13 FORMAT! 3F1C. 4, 2F 10.8) 

THE DATA HAS BEEN READ, NOW INT RUPULAT E THE TEMPERA- 
TURES AND WRITE THE LABELS. 




L= 1 

DC 410 N = 2 , 14 



LL H IU l\ - ^ f l H 

I F ( TM ( I ) .LE.TRC(N,4)) GO TO 412 





25 FORMAT (*4 IX !' THE REYNOLDS NUMBER I S ' , F 7 . 0 , / / /// ) 



ftriiciojicj o i * u \ xj i 

16 FORMAT (43X , 12 ,16X, F6 .2 ) 



1 l IX , 'TEMPERATURE ' , /) 
T SUM=0 . 0 
DC 420 J=1 , 32 
W P IT E ! 6 » 16) J , TD ( J ) 



TSLM=TSUM+TD(J) 
420 CONTINUE 
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T AVG = T SUM/32 .0 
T B ( 1 ) = T D ( 3 9 ) 

TE(20)=TD(40) 

NFITE<6,18) TAVG, TB(1 ) ,TB(20) 

18 FCRMAT (/// , 33X , 1 THE AVERAGE WALL TEMPERATURE IS 
2F6.2,' CEGREES F ',///, 35X ,' T HE ENTERING ANC EXITING 
A ' BULK TEMPERATURES ARE • ,/ ,49X , F6 . 2 , ' AND *,F6.2,///) 

CALCULATE THE LOC- MEAN TEMPERATURE AND THE OVERALL 
HEAT TRANSFER COEFFICIENT. 

C= ALOG ( ( T AVG-Tb ( I) ) / ( T AVG-T B ( 20 I ) ) 

TLM = ( TB(20)-TB(1) ) /D 

IF ( TAVG .LT .70.0) GO TC 999 

N= 1 

DC 510 K= 2 » 1 A 

IFITAVG .LE.TRCIK, 1) ) GO TO 520 
N = N + 1 

5 1 C CONTINUE 
GC TO 999 

52C E= (TAVG-TRCIN, 1 ) )/(TRC(N+l,I)-TRC(N,L)) 

RC = TRC(N,2) + E*( TRC(N+1 ,2)-TRC(N ,2) ) 
C=TRC(N,3)+E*(TRC(N+1,3)-TRC(N, 3) ) 
HAVG=(RC*C*Q*(TB(2C)-TB( 1) ) )/ (AN*TLM) 

WR ITE ( 6 , 19 ) TLM, HAVG 

19 FCRMAT (• 1* ,////, 35X, 'THE LOG MEAN TEMPERA TORE CIFFER', 
3'ENCE IS 1 ,F7.2, // ,33X, 'THE OVERALL HEAT TRANSFER*, 

5* COEFFICIENT I S • , F 8 . A , / , 3 8X , • I N BTU PER MINUTE', 

7* SCUARE FEET DEGREES F ',///// ) 

TA=(TAVG-TB( I) ) 

FM = hAVG=* : (TLM/TA) 

NR I TE (6 ,26 ) T A , HM 

26 FORMAT ( 35X ,' THE * AL L FLUID TEMPERATURE DIFFERENCE IS', 
6F6.2,// ,35X, 'THE MEAN HEAT TRANSFER COEFFICIENT IS', 

8F £ . A , // // ) 

CALCULATE THE POWER AND LOCAL HEAT TRANSFER COEFFICIENT 



TW=TAVG+459 .72 
T V = TB ( I J+459.7 2 
TTR=.00C001«( TW**4-TV**4) 

RA= .0000 1999667* AS * TTR 
DC 600 J = 1 ,16 
W=(.056£84*V(J)**2)/R(J) 

UT=U( J ) + GT (J ) 

F = U ( J ) / UT 
P(J)=W*( l.C-FJ-RA 
R 1= (209 .1518 + 1 .O/HM ) 

R2=175. 596852 
C01=A1*(T£( J)— 70.0J/RI 
CC2 = A2 * (TB(J)-7 0.0 )/R2 
CCN=CC- 1+C02 

IF(TB( J ) .LT.70.0) GO TG 999 
M = 1 

DC 610 1=2,14 

I F ( T B ( J ) .LE.TRC( I , 1) ) GO TO 620 
M = M + 1 

6 1 C CONTINUE 
C-C TO 999 

62 C G=(T8( J)-TRC(M,1) )/(TRC(M + l ,1 )-TRC(M,l) ) 
RC=TRC( M, 2 ) + G*( TRC(M+1 ,2) — TRC(M,2) ) 
C«TRC(M,3)+G*(TRC( M+l , 3 )-TRC(M, 3 ) ) 

T B ( J+ 1 ) =TB(J) + (P(J) - CCN)/ ( RC* C* C ) 



6CC 

22 



K= 2*J 



L = K-1 

TX = (TD( K ) + TD ( L ) ) /2.C 

H (J)=(RC*C*G*(TB(J+1) — T8(J) ))/(AS*(TX-T8(J) )) 

CONTINUE 

FCPMa1(//^37X, 'THE LOCAL HEAT TRANSFER CC E F F I C I E N T S ' , 
2* AND* ,// ,39X, ' THE BULK TEMPERATURES Ak E SHCWN BELCW ' 
3// / ,28X, • STATION ',5X,'DISTANCE' , 5X , ' POWER ' ,5X , 



T 
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4 'TEMPERATURE ' , 5X , ' H* , / ) 

HS=C. 0 

CC 62C N=l, 16 
HS=HS+H(N) 

SS=FLCAT ( N ) 

X(N)=-XG+SS*DX 

WRITE (6, 23) N,X(N),P!N),TB(N»,H(N) 

2 2 FCRMAT(31X,I2,7X,F8.6,4X,F8 .o,5X ,F6.2,5X,F8 .6) 

62C CONTINUE 
F A = HS / L 6 . 0 
NRITE(6,27) HA 

27 FCRMAT ( /// , 3 IX , 'THE AVERAGE LOCAL HEAT TRANSFER', 

2' COEFFICIENT IS' , F8 .4 , // / ) 

WRITE (6, 24) TB( 17) ,Tb(2C) 

24 FCRMAT(//7,41X, 'ThE FINAL BULK TEMPERATURE WAS',//, 
34 7X , ' AS CALCULATED • , F6 .2 , / / ,47 X , ' AS MEASURED ', 
4F6 .2, /////) 

GC TO 9 CO 
999 -WRITE! 6 ,29) 

29 FORMAT <// ,47X, 'SOMETHING WENT WRONG',////) 

9CC STCF 
END 
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HEAT TRANSFER C A L C U_ AT I C NS FCR A 
WALL hEATER IN A RECTANGULAR CHANNEL 



THE FLOW RATE IS 2.55 CUBIC FEET PER MINUTE 
THE REYNOLDS NUMBER IS 6C9. 



THE WALL 
STAT ION 

1 

2 
3 
A 

5 

6 
7 

e 

9 

IC 
1 1 

tl 

I A 

15 

16 

17 

18 

19 

20 
21 
22 
23 
2 A 

25 

26 

27 

28 

29 

30 

31 

32 



TEMPERATURE 

TEMPERATURE 



125.67 

125.29 
12A.71 
125.55 
125 .34 
125. A6 
12A . 29 
125.97 
12A.58 
125. 12 
12A.62 
125. A2 
125.17 
125.71 
12 A . A 5 

125. 29 

125.29 
125. 5C 
12A. 92 
125. OC 
125. A6 
12 A. 87 
125.21 
125.38 
12A. 96 
125. OA 
12A. A5 
125. C£ 
12A . A 1 
12A . A5 
125. 12 



THE AVERAGE WALL TEMPERATURE IS 125. C8 DEGREES F 

THE ENTERING AND EXITING EULK TEFFERATURES ARE 

80. 13 AND 90.87 



148 



35.34 



THE LOG MEAN TEMPERATURE DIFFERENCE IS 



THE OVERALL HEAT TRANSFER COEFFICIENT IS 

IN 6TU PER MINUTE SQUARE FEET CEGREES 



0.C460 

F 



THE WALL FLUID TEMPERATURE DIFFERENCE IS 44.95 
THE MEAN HEAT TRANSFER COEFFICIENT IS 0.04C3 



THE LOCAL FEAT TRANSFER COEFFICIENTS ARC 
THE BULK TEMPERATURES ARE SHCWN BELOW 



T ION 


DISTANCE 


POWER 


TEMPERATURE 


H 


1 


0.011416 


0.088965 


8C . 13 


0.132104 


2 


0.034249 


0.038455 


£2 . C8 


0.C59C25 


3 


0.057082 


0 .043 506 


8 2.52 


C .C67486 


4 


0. 079915 


0.C33487 


83 .86 


0.053317 


5 


0. 1C2748 
0. 125581 


0 .028635 


84.58 


0.045572 


6 


0 .026755 


85.20 


0.C43511 


7 


0. 14S414 


0.027666 


£ 5.77 


C .045427 


8 


0.171247 


0.021361 


86.36 


C . C34521 


9 


0. 194030 


0.018317 


86.81 


0.029719 


10 


0.216913 


0.020845 


87. IS 


C .034529 


11 


0.239746 


0 .0 ic224 


8 7.62 


0 .C26526 


12 


0. 262575 


0 .020555 


87.55 


0 .034747 


13 


0.2854 12 


0 .0 15C3 1 


£8.38 


C.C248C7 


14 


0.308245 


0 .013437 


66 .69 


C .C22268 


15 


0.331C78 


0.G1 1453 


88.55 


0.018718 


16 


0.353911 


0 .019532 


£5 . 18 


0.C34C46 



THE AVERAGE LCCAL HEAT TRANSFER COEFFICIENT IS 0.C442 



THE FINAL BULK TEMPERATURE WAS 
AS CALCULATED 89.58 
AS MEASUREC 90.87 
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HEAT TRANSFER CALCULATIONS FCR A 
WALL HEATER IN A RECTANGULAR C F AN N EL 



THE FLGW RATE IS 3.66 CUEIC FEET PER MINUTE 



THE REYNOLDS NUMBER IS 874 . 



THE WALL TEMPERATURE 



STATION 


TEMPERATURE 


1 


131.99 


2 


133.24 


a 


132.86 


4 


132.07 


5 


133 .32 


6 


132.82 


7 


133.94 


8 


131.45 


9 


134. 52 


10 


131.99 


11 


133.24 


12 


132. 74 


13 


133. 82 


14 


133 . 15 


15 


133.82 


16 


132.37 


17 


132 .99 


18 


133. 2£ 


19 


133.24 


20 


132.45 


21 


132.03 


22 


132. 53 


23 


132 .4 1 


24 


132.74 


25 


133.49 


26 


132 .78 


27 


133.15 


28 


132.4 1 


29 


133.36 


30 


132.24 


'31 


132.61 


32 


133.03 


THE AVERAGE WALL 


TEMPERATURE IS 132.88 DEGREES F 


THE ENTERING AND EXITING EULK TEMPERATURES ARE 




76.99 AND 88.51 
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49. SI 



THE LOG MEAN TEMPERATURE DIFFERENCE IS 

THE OVERALL HEAT TRANSFER COEFFICIENT IS 

INI BTU PER MINUTE SQUARE FEET CEGREES 



0 .0551 
F 



THE WALL FLUID TEMPERATURE DIFFERENCE IS 55.89 
THE MEAN FEAT TRANSFER COEFFICIENT IS 0.C452 



THE LOCAL FEAT TRANSFER COEFFICIENTS AND 
THE BULK TEMPERATURES AkE SFCWN BELOW 



STATION 


DISTANCE 


POWER 


TEMPERATURE 


H 


I 


0.011416 


0 . 123292 


76.99 


0. 149215 


2 


0.034249 


0 .063262 


78.68 


0.C78692 


3 


0. 057C82 


0.064357 


79.64 


C .080715 


4 


0.079915 


0 .056222 


£ C . £ 2 


C. C72119 


5 


0.102748 


0.045533 


8 1.67 


C . C584C3 


6 


0. 12556 1 


0 .046022 


£2.26 


C .060099 


7 


0. 148414 


0 .046040 


82 .C6 


C .060267 


8 


0. 17 1247 


0.C36775 


83.76 


0 .0438 15 


9 


0. 194060 


0.C316C5 


34.21 


C.C42C99 


10 


0.216913 


0 .036077 


84 .79 


C.C49CC4 


1 l 


0.239746 


0.025993 


85.33 


0.035561 


12 


0.262579 


0.036541 


85.71 


C .050813 


13 


0.285412 


0 .023612 


86.26 


0 . C39289 


14 


0. 30 £245 


0.027 137 


66.69 


0.037734 


15 


0.331078 


0 .022255 


6 7.C9 


0.030922 


16 


0.35391 1 


0 .029529 


87.42 


0 . C41776 


THE 


AVERAGE LOCAL 


HEAT TRANSFER 


COEFF IC IENT 


IS 0 .0585 



THE FINAL BULK TEMPERATURE WAS 
AS CALCULATED 87.85 
AS MEASUREC 38. 51 
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HEAT TRANSFER C A L CU L AT I Ci\ S FCR A 
WALL HEATER IN A RECTANGULAR CHANNEL 



THE FLOW RATE IS 4.88 CUBIC FEET FER MINUTE 
THE REVNGLDS NUMEER IS 1166. 



THE WALL 
STATION 

1 

2 

3 

4 

E 

6 

7 

8 

5 

10 
1 1 
12 

13 

14 

15 

16 

17 

18 
19 
2 C 
2 1 
22 

23 

24 

25 

26 

27 

28 
25 
3C 

•31 

32 



T EMFERATURE 

TEMPERATURE 

143 . 16 
144. 8 C 

144.06 
143.52 
144. 22 

143. 6 1 
144.35 
142.34 
144.63 

142.38 
143.81 

143.07 
143.73 

143.16 
143. 4C 
141.80 

142 .38 
142.87 

143.07 

142.38 
142.54 
143 .28 
142 . 46 
14 ^ . 42 
143.52 
142. 7C 
143 . £8 
142.66 

143. 16 
142. 05 
142 .42 
142.87 



THE AVERAGE WALL TEMPERATURE IS 143.13 DEGREES F 

THE ENTERING AND EXITING BULK T E N F ERATUR ES ARE 

79.82 ANC 51.96 
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THE LOG MEAN TEMPERATURE DIFFERENCE IS 57.02 

THE OVERALL HEAT TRANSFER CCEFFI Cl EM IS 0.0666 
IN BTU PER MINUTE SQUARE FEET CEGREES F 



THE WALL FLUID TEMPERATURE DIFFERENCE IS 63. 3L 
THE MEAN FEAT TRANSFER COEFFICIENT IS O.C6CO 



THE LOCAL FEAT TRANSFER COEFFICIENTS ANC 
THE BULK TEMPERATURES ARE SHCWN BELCW 



TI CN 


DI STANCE 


PCWER 


TEMPERATURE 


H 


1 


0 . 0 1 1 4 1 6 


0 . 14793C 


79 .62 


0. 155C66 


2 


0. 034249 


C.C81882 


81 .53 


0 .087870 


3 


0.057C82 


0.080664 


62.47 


C. C87608 


4 


0.079915 


0 .070126 


83 .39 


0 .077752 


5 


0. 1C2748 


0.054745 


84.15 


0.061034 


6 


0. 125581 


0 .059 155 


64 . 62 


C.C66625 


7 


0.148414 


0 .C55656 


85.50 


C .063448 


8 


0. 17 1247 


0 . C44 4 5 1 


86. 13 


C.C5157L 


9 


0 . 194080 


0 .037650 


86 .64 


C. 04 3745 


1C 


0.216913 


0.044400 


87 . C6 


0.052157 


11 


0.235746 


0 .04023 1 


87.56 


0 . C47309 


12 


0.262579 


0 .045 171 


86 .02 


C . C542C4 


13 


0.285412 


0 . 0340 27 


88.53 


0.040176 


14 


0.306245 


0 .038087 


86.51 


C .045604 


15 


0.331078 


0 .026458 


85 .34 


C .031473 


16 


0.353511 


C.C35C7C 


65.63 


0.C42578 


THE 


AVERAGE LCCA-L 


HEAT TRANSFER 


COEFFICIENT 


IS 0.063 



THE FINAL BULK TEMPERATURE WAS 
AS CALCULATEC 90.03 
AS MEASURED 91.56 
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HEAT TRANSFER CALCULATIONS FCR A 
WALL HEATER IN A RECTANGULAR CHANNEL 

THE FLOW RATE IS 6.10 CUBIC FEET FER MINUTE 
THE REYNOLDS NUMBER IS 1457 . 



THE WALL TEMPERATURE 
STATION TEMPERATURE 



1 

2 

3 

4 

5 

6 
7 

e 

s 

i? 

12 

13 

14 

1 5 

16 
17 

ie 

19 

2 C 
21 
22 

23 

24 

25 

26 

27 

28 
25 

30 

31 

32 



71 

c; > 



139. 

140. 
139. 

138 

139 
139. 

140 
137 
140 
137 
139.0 
138.26 
139 



71 

8 C 
42 
2 5 

4 1 
51 
61 

■ 8 * 



138. 
135 , 

137 
138, 
138, 

139 

137, 

138. 

138 
138 
138. 
L39. 
138 
139.34 
138.42 

139.25 
137 . ee 

138.25 
139,05 



13 
17 
21 
47 
05 
71 
0 5 
88 
13 
88 
22 
2 
4 
51 



THE AVERAGE WALL TEMPERATURE IS 138.84 DEGREES F 

THE ENTERING AND EXITING BULK TEMPERATURES ARE 

81.23 AND ST. 57 
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THE LOG MEAN TEMPERATURE CIFFERENCE IS 52.27 

THE OVERALL HEAT TRANSFER CCEFFICIENT IS 0.C779 
IN BTU PER MINUTE SCLARE FEET CEGREES F 



JHE WALL FLUID TEMPERATURE DIFFERENCE IS 57. 6 L 
ThE MEAN HEAT TRANSFER CCEFFICIENT IS O.C7C7 



THE LOCAL FEAT TRANSFER COEFFICIENTS AND 
THE BULK TEMPERATURES ARE SHCWN E E L C W 



T I CN 


DISTANCE 


POWER 


TEMPERATURE 


H 


I 


0.011416 


0 .148667 


8 I .23 


0. 169629 


2 


0.034249 


0.C82567 


82 .60 


0 .C97266 


3 


0 . 05 7082 


0.082981 


82.26 


C.C98860 


4 


0.079915 


0 .073819 


84 . 12 


0.C89528 


5 


0. 102748 


0.C59164 


84 . 80 


0.071915 


6 


0 . 125581 


0.060822 


85.24 


C .C75466 


7 


0.148414 


0 .Gr>8258 


85.50 


C .C72912 


8 


0. 17 1247 


0. C496C2 


86.43 


C .C62756 


9 


0. 1 940 80 


0 .0429 77 


86 .88 


C .C544SC 


10 


0.216913 


0.049149 


87.27 


0 .062939 


u 


0.239746 


0.043822 


cl . i2 


0.C562C6 


12 


0.262579 


0 .049fc7 1 


88.11 


C .C645 17 


13 


0.285412 


0 . 039 2 6 C 


88.56 


0 .050447 


14 


0. 308245 


0 .044015 


88.52 


C.C57253 


15 


0.331078 


0 .031205 


89.22 


C .C40521 


16 


0.35391 1 


0.0395(6 


£ 5 . 6 C 


0.052026 


THE 


AVERAGE LOCAL 


HEAT TRANSFER 


COEFFICIENT 


IS 0.C73- 



THE FINAL BULK TEMPERATURE NAS 
AS CALCULATEC 89.95 
AS MEASURED 91.57 
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HEAT TRANSFER CALCULATIONS FCR A 
WALL FEATER IN A RECTANGULAR CHANNEL 



THE FLOW RATE IS 6.10 CUBIC FEET FER MINLTE 



THE REYNOLDS NUMBER IS 1457. 



THE WALL TERFERATURE 



AT ION 


TEMPERATURE 


1 


132. 66 


2 


133.52 


3 


132.37 


4 


131.29 


5 


132. 12 


6 


131.54 


7 


132 . 49 


8 


129.79 


9 


132.66 


1C 


130.04 


l 1 


131.37 


12 


130.62 


13 


131.37 


14 


130. 54 


15 


131 .4 l 


16 


129.71 


17 


130.4 1 


18 


130.95 


19 


13i .24 


20 


120.46 


21 


130. 5C 


22 


131.16 


22 


130. 4 1 


24 


130.54 


25 


131.83 


26 


131. OE 


27 


131.83 


28 


131 .00 


29 


131. 87 


30 


130.56 


31 


131 . 12 


32 


131. 7C 



THE AVERAGE WALL TEMPERATURE IS 131.25 DEGREES F 



THE ENTERING AND EXITING SULK TEHFERATURES ARE 

78.32 AND 88.60 
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THE LOG MEAN TEMPERATURE LIFFERENCE IS 47.61 

THE OVERALL HEAT TRANSFER COEFFICIENT IS 0.C861 
IN BTU PER MINUTE SCUARE FEET CEGREES F 



THE WALL FLU 10 TEMPERATURE DIFFERENCE IS 52.94 
THE MEAN HEAT TRANSFER COEFFICIENT IS 0.0775 



THE LOCAL FEAT TRANSFER COEFFICIENTS ANC 
THE BULK TEMPERATURES ARE SHCWN BELCW 



T ION 


DISTANCE 


POWER 


TEMPERATURE 


H 


1 


0.011416 


0.148165 


78.32 


0. 182080 


2 


0 . 034249 


0. C8 175C 


79.68 


C . 104942 


3 


0 .057082 


0 .082228 


8 C .43 


C . 107013 


4 


0. C79915 


0.C73501 


81.18 


0 .098191 


5 


0 . 102748 


0.060481 


8 1.86 


0.081235 


6 


0.125581 


0 .061153 


82.41 


C .C836 17 


7 


0. 148414 


0. C58306 


82.97 


0.080580 


8 


0. 171247 


0 . 050CCC 


83 . 5C 


C. 070163 


9 


0.194080 


0 .044486 


83.96 


0.C62616 


LO 


0.216912 


0.051C 15 


8 4.36 


0 .072384 


11 


0.239746 


0 .04^598 


84 . 82 


G . C63653 


12 


0.262579 


0 .050545 


85.23 


0 .C73555 


13 


0.285412 


0.C4UC9S 


85.69 


C.C57226 


14 


0.306245 


0 .045 ice 


86 . C 5 


0 .C6516C 


15 


0.331078 


0.C34332 


86 .46 


0 . C49 68 6 


16 


0.3539 1 1 


0 .035976 


8 6.77 


C. 052272 


THE 


AVERAGE LOCAL 


HEAT TRANSFER 


COEFFICIENT 


IS 0.C815 



THE FINAL BULK TEMPERATURE WAS 
AS CALCULATEC 87. OS 
AS MEASURED 88. 6C 
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